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NOTE: 


Generally, a variety of watershed data are compiled on a 
calendar year basis. However, the water year, beginning 
October 1 and ending September 30, has proven best for 
hydrologic comparisons. 
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INTRODUCTION 


Cooperative watershed research between the Agricultural Research 
Service, U. S. Department of Agriculture, and the Bureau of Land 
Management, U. S. Department of Interior, was initiated in 1968 
under Cooperative Agreement No. 14-11-0001-4162(N). Also, the 
Memorandum of Understanding, dated July 6, 1960, which is part of 
the Cooperative Agreement, specifies the overall responsibility of 
each agency. 


This interim report summarizes progress and results on the Reynolds 
Creek Watershed and supporting studies on the Boise Front from 
October 1 through September 30, 1981. Data collection, processing, 
analyses, and reporting are according to the FY 1981 work plan. 
Progress reports are given by the individual sections of the work 


plan. A copy of the FY 1981 work plan precedes the progress 
reports. 


Supporting information and data are presented in Northwest 
Watershed Research Center Annual Reports for 1972 and prior years, 
anc I VertneertmanevOr Ge NGOs 1," yuo. 4s 5 Oe le esos Ose atid 
for the ARS-BLM studies in the Reynolds Creek Watershed under 
Cooperative Agreement No. 12-14-5001-6028. 


The research program at the Northwest Watershed Research Center, 
which includes field, laboratory, and computer-based studies, 
services action agency needs and contributes to the sciences of 
watershed hydrology and watershed management. Cooperative research 
With user agencies, such as BLM and SCS, and with State 
Universities in Idaho and Utah serves to enhance the application of 
research findings. 


Research investigations at the Center include the following: (not 
in order of priority) 


1. Snow Hydrology - Investigations on deposition and 
redistribution of snow, energy balance and melt processes of 
Snowpacks, and water supply forecasting of snowmelt runoff. 
Research on improving water supply forecasting utilizes SCS 
SNOTEL and snow course data, and is conducted in cooperation 
With them and the National Weather Service. Snow research also 
contributes to the ARS AgRISTARS program. Snow research for 
mid—-elevation rangeland watersheds is conducted on the Reynolds 
Creek Experimental Watershed. (A. L. Huber and K. R. Cooley, 
Hydrologists). 


2. Water Quality - Investigations center on the effects of grazing 
systems on rangeland water quality. This research has 
identified rangeland non-point pollution sources and is 
evaluating the effects of management practices. Winter holding 
and feeding areas in the Reynolds Creek valley are being 
studied as a potential pollution source. The occurrence and 
movement of fecal coliform bacteria are being studied on the 
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Reynolds Creek Watershed and also on the Boise Front. Computer 
modeling of streamflow water quality is underway. Cooperative 
research is with BLM, the University of Idaho, and the Idaho 
Fish and Game Department. (G. R. Stephenson, Geologist). 


Runoff and Water Yield - Investigations are developing and 


evaluating watershed runoff and water yield models. Winter 
season runoff events are major flood and sediment producing 
events on range and agricultural watersheds of the interior 
Northwest. Rainfall, occurring with melting snow on frozen 
soil, combines to generate these critical runoff conditions. 
Research on frozen soil hydrology is cooperative with the 
University of Idaho and other Pacific Northwest Area locations. 
(C. W. Johnson, Research Hydraulic Engineer, J. P. Smith, 
Hydrologist, C. L. Hanson, Agricultural Engineer, and K. R. 
Cooley, Hydrologist). 


Erosion and Sediment Yield -— The Universal Soil Loss Equation 


(USLE) requires modifications when applied to semiarid 
rangeland watersheds in the Northwest, especially for use with 
frozen soil and melting snow conditions. The USLE and sediment 
yield equations, such as the MUSLE, used by BLM and SCS on 
rangelands are being evaluated. Bedload studies are 
identifying the magnitude of bedload and suspended sediment 
transport in steep, rocky channels. Streams on the Reynolds 
Creek Experimental Watersheds are utilized in these erosion and 


sediment studies. (C. W. Johnson, Research Hydraulic 
Engineer). 


Rangeland Vegetation - Investigations on the forecasting of 
rangeland forage production are underway. In cooperation with 
BLM, vegetation study sites established on the Reynolds Creek 
area are being used to evaluate the effects of livestock 
grazing. Also, rangeland forage forecasts were developed in 
cooperation with BLM in Montana. On the Reynolds Creek 
Watershed, three plant adaptability nurseries were established 
in cooperation with the Forest Service. A coordinated ARS 
project on rangeland resource modeling is underway. (J. R. 
Wight, Range Scientist, and C. L. Hanson, Agricultural 
Engineer). 


Precipitation -— Research is focused on developing and testing 
models for predicting the precipitation input to rangeland as 
influenced by elevation, aspect, and season. A study on "acid" 
precipitation is being initiated. (C. L. Hanson, Agricultural 
Engineer, and G. R. Stephenson, Geologist). 











7. Soils - Research on modeling soil infiltration is being 
conducted in cooperation with SCS. Published soil water 
retention data have been used to predict Green—Ampt 
infiltration equation parameters. Procedures to predict the 
influence of management practices on infiltration and soil 
Water storage are under development. (D. L. Brakensiek, 
Research Hydraulic Engineer, and W. J. Rawls, Hydrologist - 
Beltsville, Maryland). 





8. Watershed Modeling - Combining all components of rangeland 
watershed hydrology and resources into a rangeland watershed 
model is underway. Development of a watershed model will 
provide the manager with a tool for making rangeland management 
decisions. Alternative watershed practices and conditions can 
be evaluated without lengthy and costly field installations. 
Component models under development and/or testing include 
watershed evapotranspiration, soil moisture Storage, and 
infiltration and runoff. A major effort underway involves the 
coordination of major ARS activities in Rangeland Hydrology and 
Resource Modeling: (C.°L. Hanson, Agricultural Engineer, J. R. 
Wight, Range Scientist, and Staff). 


Introduction Figure 1 locates major experimental sites on Reynolds 
Creek Watershed. 


Additional copies of this report or further information on reported 
work can be obtained from: 


Northwest Watershed Research Center 
USDA, Agricultural Research Service 
1175 South Orchard, Suite’ 116 
Boise, Idaho 83705 


REYNOLDS CREEK Recet DATA NETWORK 
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Figure 1.--Reynolds Creek Watershed. 
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BLM-ARS ANNUAL WORK PLAN FOR FY 1981 


INTRODUCTION: The following work plan items contribute to the 
objectives of Paragraph III of the Bureau of Land Management 
Interagency Agreement No. YA-515-IA8-21 dated September Too Code 
Certain items in the work plan represent a continuation and/or 
completion of work from previous years. In some cases, the 
continuation represents collection of hydrologic data or watershed 
resource inventories. This is required for sampling the influence 
of climatic variability on hydrologic factors and cumulative 
effects of rangeland management practices. USDA-ARS watershed 
management research not a part of this work plan will, in many 


instances, supplement and/or complement the following work plan 
items: 


TegetPRECIPLTATION 


Continue operation of the precipitation network on Reynolds 
Creek. A new site (098X97) was established at the Nancy runoff 
plot on October 1, 1980 and the old Summit site (049X61) was 
reactivated on the same date. The four sites on the Boise 
Front will be discontinued on January 1, 1981. Monthly and 
annual precipitation models will be applied to comparable 
rangeland areas in Idaho, Oregon, and Nevada. Stochastic 


modeling of precipitation events on Reynolds Creek will 
continue. 


ce» VEGETATION 


Complete testing of Saxton soil water model. Initiate an 
analysis of Reynolds Creek data on species composition and 
cover at nine grazed and four ungrazed sites to determine 
Sampling variability with application to the design of 
vegetation monitoring systems. Bitterbrush utilization surveys 
will continue on the Boise Front, depending on technical 
Support from Idaho Fish and Game. Finalize the preparation of 
a Reynolds Creek Watershed soil/vegetation data base. 


3. RUNOEP 


Continue runoff data collection at two main stem, two tributary 
and three source watersheds on Reynolds Creek. Runoff sites on 
the Boise Front study areas will be terminated January 1, 1981. 
Runoff records from the Saval Ranch (Nevada) project will be 
processed as funding and records are available. Runoff model 
testing will commence with Reynolds Creek Watershed data. 

Green and Ampt infiltration parameters will be utilized in 


modifying runoff prediction equations, such as the SCS runoff 
equation. 


EROSION AND SEDIMENT 


Continue data collection at four Reynolds Creek sites and 
terminate the two sites on the Boise Front on January 1, 1981. 
An evaluation of the relationship between plant biomass and 
vegetative cover on ungrazed sites will be conducted to improve 
the cover factor in USLE. Continue the evaluation of USLE, 
MUSLE, and PSIAC factors for rangeland watersheds. 


WATER QUALITY 


Continue the evaluation of water quality models with Reynolds 
Creek data. Determine the influence of streamflow channel 
characteristics, such as sediment, temperature, and flow 
velocity on quality indicators, such as fecal coliform. 
Continue water quality observations on the Boise Front 
rest-rotation system and compare the effects of cattle and deer 
grazing. Investigate the survival of E. coli in stream water, 
in stream sediments, and in off-stream sites. 














SUMMARY OF WORK PLAN RESEARCH FINDINGS 


BY WORK PLAN SECTIONS 
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Responsible for data reduction and 


processing. 


Assists in hydrologic data 
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Reynolds Creek precipitation: The precipitation network in operation 
on the Reynolds Creek Watershed is shown in Introduction, Figure 1. 


The Tour precipitation sites: Listed MntTable 121 represent the 
precipitation conditions that existed on the watershed in water year 
1981. Precipitation ranged from 2.3 inches below average at the low 
elevation site 076X59 (3965 feet) to 14.5 inches below average at the 
high elevation site 176X07 (6760 feet). The much below average 
annual precipitation amounts were due to the very dry conditions 
during the fall and early winter. Precipitation during the spring 
and early summer (March-June) was at or above average, which resulted 
in good vegetative growth. The annual total at the Boise airport was 
0.7 inch above average because of above average precipitation during 
March and April and near normal precipitation most other months. 


Boise Front precipitation -— 1981 


Precipitation gages were removed from the Boise Front the week 
following January 1, 1981. 


Wyoming shield gage study: Preliminary data analysis indicates that 


rain and wet snow catch by the Wyoming shield gages is equivalent to 
the catch by the dual-gage system when 1.8 is used as the coefficient 
in the dual-gage amount equation. We are continuing this study 
because our analysis of snowfall measurement has not been completed. 


Heated precipitation gages: A study was initiated during the fall of 
1980 to determine if the snow catch by heated precipitation gages is 
comparable with that measured by standard weighing and recording 
gages. Preliminary analysis indicates that the heated gages catch 
about 30 percent less. This study is continuing during the 1981-82 
winter because of the small number of snow events during 1980-81. 


Saval Ranch: ARS involvement in the hydrologic aspects of the 
BLM-ARS Saval Ranch Project continued this year. Major emphasis 
centered on insuring that quality hydrologic data would be available 
when model development and testing is initiated. To meet this 
objective, efforts were made to improve and standardize precipitation 
collection equipment, a data processing technician was hired and is 
being trained, and a computerized data file was initiated. 


During the week of June 22-26, 1981, Keith Cooley and Dave Robertson 
visited the Saval Ranch for the purpose of installing better gage 
mounting brackets at all of the precipitation sites, and to tighten 
and adjust the Wyoming Shields as needed. The original wooden 
mounting brackets had become loose due to drying and splitting caused 
by weathering processes and vibration. The new steel brackets 
consisted of 1/4" thick steel plates bolted directly to the bottom of 
the gage (eliminating the 3-legged stand provided by the 
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Table 1.1.--Water year precipitation (inches) at four locations on Reynolds Creek Watershed and the Boise At poKES as 
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Site Se Year Oct Nov Dec Jan Feb Mar Apr May June July Aug Sep Total 
ft) 2. 
076X59 3965 1981 -350 +884 1.528 - 283 888 1.090 1.722 1.021 .930 .020 .070 .090 8.876 
1963-1981 842 1.176 L203 1.493 . 787 .917 oi eH - 885 1.338 Sate} .746 .569 11.190 
116X91 4760 1981 -707 1.074 2.370 -640 1.943 1.574 2.617 1,599 1.445 .130 .110 .190 14.399 
1963-1981 1.434 2.054 2.341 2.666 1.545 1.682 1.775 1.260 1.601 389 AOS .796 18.246 
155X07 5410 1981 died 1.586 3.864 1.307 2.794 3.432 2.891 2.228 1.226 .290 .120 3 fei) 20.914 
1963-1981 1.985 3.415 3.746 4.552 2.700 2.784 2.357 1.771 1.836 . 608 1.034 1.029 Phos Wi 
176X07 6760 1981 1.083 ke ft) 4.061 2.830 2.820 5.124 3.039 2.767 1.457 -310 -150 ~458 26.376 
1963-81 2.219 5.365 5.843 7.880 4.396 4.264 3.514 2.429 2.180 seit) P1020 ae lg 40.886 
Boise 

Airport 2838 1981 . 300 1.260 1.490 1.200 1.020 2.760 1.930 .950 .770 . 230 .130 . 360 12.400 
1941-1981 .829 1.316 1.351 1.480 1.143 1.118 1.167 1.233 1.000 . 206 7339 pepe 11.695 
ee see ee ee ee ee oe, ee ee OS Ce 
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Raingage locations are shown on Introduction, Figure 1. 





manufacturer), three 1/2" diameter leveling bolts, and three angle 
braces anchored to the posts with lag bolts. Overall, the new 
brackets appear to have increased mounting sturdiness, and should 
reduce gage vibration caused by the wind. 


The gage orifices were also repositioned to be in the same plane as 
the top of the inner shield of the double shielded Wyoming Shield 
configuration. Some were found to be up to 6 inches low, or high, 
probably due in part to settling of the shields since installation. 
Although gage catch may not have been affected Significantly by the 


discrepancy in height, it seemed best for future comparisons to have 
them all positioned the same. 


The calibration of all gages was checked, and all appeared to be 


working satisfactorily, although some were off Slightly at the higher 
depth readings. 


In October, Randy Evans was hired to process the Saval data. He has 
been trained to mark the precipitation charts, and has marked all of 
the charts from June, 1979 to September, 1981. He is now in the 
process of being trained to digitize the chart data, and enter it 
into the computerized data file that is being developed. 


Neil Hutten, stationed at Elko, Nevada, also received training in 
precipitation gage calibration and field procedures during a two-day 
visit to Boise in late August. 


The computerized data file that is being developed will contain all 
of the hydrologic data being collected at Saval Ranch. This will 
include breakpoint continuous precipitation and runoff from the small 
watersheds; and daily values of evaporation, maximum and minimum 
temperature, solar radiation, humidity, and runoff from the four USGS 
gaging stations as processed and made available by USGS. The 
precipitation, small watershed runoff, and evaporation charts will be 
processed at Boise. Daily solar radiation, maximum and minimum 
temperature, and humidity will be recorded on data sheets by Neil 
Hutten and sent to Boise, so these data can be entered into the file. 
Ouput from the data file can then be obtained in several formats, 
including reproductions of the file on magnetic tape, and analyzed 
information such as storm EI, seasonal or annual EL, intensity 
characteristics, etc. 


Examples of breakpoint precipitation data for three storms were 
presented in Interim Report No. 11. An example of daily and monthly 
precipitation data is presented in Table 1.2 for the last half of 
1979 for station #4. Table 1.2 represents a standard type of output 
that will be available, but almost any format could be provided. 


During the week of November, 9-13, Keith Cooley participated in a 
modeling workshop at Elko to help develop an initial version of an 
integrated model for the Saval Project. Many of the hydrologic 
components were based on previous ARS research, and were described by 
Cooley for inclusion in this model. 


ae 


At both the June and November visits to Elko, Keith Cooley 
participated in sessions held to discuss future direction of the 
hydrologic aspects of the Saval Project. Recommendations were made, 
based on previous watershed hydrology experience at Boise, to support 
and facilitate the collection of the proper hydrologie data that will 
be required to meet future project objectives. 


Table 1.2.--Daily and monthly summary of precipitation for July through 
December, 1979; station #4, Saval Ranch project, Elko, Nevada. 


MONTH 
DAY JUL AUG SEP OCT NOV DEC 
ah -——— —-——— -_—— ——— -—— ——— 
2 aoe — — = mee 0.020 
3 ae a a ~~ 0.290 --- 
4 —— ae ees —— 0.120 —— 
5 ae = oe ee sey eae 
6 mee ae ee peers = —— 
a aia an ass oes a aa 
8 Jagat oe a SEIS Aes ae 
9 aed —— — =< om Ree, 
10 _— --- --- --- --- --- 
11 --- --- --- --- --- =e 
12 --- 0.039 --- --- --- --- 
13 i. 0.091 --- -—— --- — 
14 --- 0.220 --- --- --- --- 
15 -—- —- --- 0.140 ~-- $22 
16 -—- 0.010 --- —- ~-- = 
17 —— 0.130 --- --- 0.100 2s 
18 --- 0.160 --- 0.151 0.120 — 
19 --- 0.090 --- 1.697 0.010 — 
20 --- 0.020 --- 0.282 a ox 
21 0.300 0.010 --- 0.130 --- 0.140 
22 0.700 --- --- aes 0.226 0.040 
23 0.010 --- -—~ --- 0.154 0.100 
24 --- --- --- --- 0.330 see 
25 --- --- 0.390 --- 0.590 0.100 
26 --- --- --- — 0.050 ai 
27 --- -—— -~- es pe ee 
28 —= --- = 0.080 a a 
29 --- --- --- 0.070 --- ot 
30 -—- 0.150 ~-- 0.055 = — 
31 --- --- —-~ 0.135 --- cree 
TOTAL 1.010 0.920 0.390 2.740 1.990 0.400 
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ANNUAL AND MONTHLY PRECIPITATION GENERATION 
FOR IDAHO AND SURROUNDING STATES 


Data from 33 stations in Idaho and 25 stations in states surrounding 
Idaho are being used to develop procedures for generating annual and 
monthly precipitation amounts. Most station data sets being used are 


4O year (1940-79) records; however, a few records are shorter because 
of missing data. 


The procedure that was developed to generate annual series is based 


on the lognormal distribution and is a shortened version of the 
procedure discussed in Interim Report No. 11. Mean annual 


precipitation is the only variable required in this procedure because 
only mean annual precipitation and the associated coefficient of 
variation are used in the generation. 

The equation for generating annual precipitation using the lognormal 
G2iUstributions 1s: 


sO exp(Up + pre [1] 


where, 
Y = generated annual precipitation values 
py = the mean of the log, values 


Ofny = the standard deviation of the log, values of the site 
where data are available 


y = a pseudo-random normal deviate [N(0,1)]. 


The equations for obtaining u 


enx and ° pny from the mean annual are 
(Haan, 1977): 


— wp ary 2 2 
Uprny = LenlX /(C, + 1)] [2] 
x 2 5 
Oar Cenc, + 1)] [3] 
and, 
Cre 0, 262 07 00018% [4] 
(R2= 0.613) 
where, 
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X = average annual precipitation (mm) obtained from records or maps 


sa = eoefficient of variation. 


The generating procedure is as follows: 


1. Compute the mean annual precipitation (X) for the location in 
question from climatic records, maps, etc. 


2. Compute the coefficient of variation by using Equation [4]. 
3. Compute Up, and Spy, by using Equations [2] and [3]. 


4, Obtain as many values of y as is required from a table of 
pseudo-random normal deviates [N(0,1)] or from computer programs for 
generating these values. 


5. Compute the annual precipitation series by using Equation [1]. 


Discussion: The coefficient of variation (C_) values were plotted on 
maps, to determine if it varied by location Yn the study area. No 
regional trends could be found from the maps, so all of the data were 
analyzed in one group. There is considerable scatter among the data 
points, but the Coefficient of Determination (R72 = 0.613) indicates 
that there is a reasonable relationship between average annual amount 
and the Coefficient of Variation. . 


The summary of 50-year simulations listed in Table 1.3 are for 
locations where the average annual precipitation was calculated to be 
approximately 250 and 750 mm. The simulated standard deviation was 
greater than measured by about 12 percent. The simulated range of 
annual totals was about 12 percent greater than the measured values. 


MONTHLY PRECIPITATION GENERATION 


The monthly generation procedures currently being developed for Idaho 
and surrounding states are based on the monthly precipitation scheme 
developed for the Reynolds Creek Watershed data set, aS reported in 
Interim Report No. 11. Preliminary analysis indicates that the 
generating procedures outlined in Interim Report No. 11 can be used 


to generate monthly precipitation amounts for Idaho and surrounding 
States, 


REFERENCES 


Haan, C.T. 1977. Statistical methods in hydrology. The Iowa State 
University Press, Ames, Iowa. 378 pp. 
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Table 1.3.--A summary of 50-year simulations where the calculated average 
annual precipitation used in equations [1] through [4] was 
approximately 250 and 750 mn. 








Average Coefficient of Standard Range of 
Precipitation Variation Deviation Annual Totals 
(mm) (mm) (mm) 

Simulated 251 Wee 61 157 =» 425 
Deer Flat Dam, ID 251 wee 53 126° = e335 
Malheur, OR 247 eae s ss ta 131 “s° 3356 
Simulated 747 18 137 476 - 1162 
Saint Maries, ID 762 S/ 126 465 - 1021 


Kellogg, ID 773 .16 125 485 - 1037 
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Reynolds Creek (Reynolds Creek site locations are shown in 
Introduction, Figure 1). 


Soil water balance models: Progress made in testing Wight's and 
Saxton's soil water balance models consisted of: (1) completion of a 
test data file containing six or more years of data for each of four 
Reynolds Creek watersheds; (2) transforming the format of Saxton's 
model to make it compatible with the Center's computer, and checking 
it against a test data set; (3) obtaining the necessary input factors 
and boundary condition parameters to make a calibration run with 
Saxton's model using the 1979 data from the Flats Site; (4) running 
Wight's model continuously for six years (1975-1980) using the 
calibration year; and (5) making a comparison of soil moisture output 
from the models with measured soil moisture at the Flats site for 


WM ee 


The computerized hydrologic data file for testing the soil water 
balance models has been completed. The file consists of daily values 
of precipitation, runoff, solar radiation, and maximum and minimum 
temperature for six or more years at four Reynolds Creek sites. 

Daily values of solar radiation for the Boise airport for the same 
years have been obtained and entered into the file for correlation 
and use in estimating missing data. Daily pan evaporation data for 
the summer periods at the two upper watersheds (Reynolds Mountain and 
Lower Sheep Creek), and at a third site, representing the two lower 
watersheds (Nancy and Flats), is also included. Soil moisture data 
for five layers, down to about 136 em, are used as a check on the 
models. The four watersheds vary in elevation from 1180 m to 2160 mM, 
and in size from 0.8 ha to 40 ha. The annual precipitation ranges 
from 24 cm at the lower sites to 102 em at the upper site. 

Vegetation consists mainly of grasses, Sagebrush, and shadscale at 
the Flats, and grasses and sagebrush at the other Sites. However, 
they do represent an increasing amount of larger shrub and woods with 
increasing elevation. Soils vary from fine loam to gravelly loam. 
These four sites were chosen for testing the models because of data 
availability and the range site characteristics represented. 


Although the copy of Keith Saxton's model provided to us was written 
in FORTRAN, it had to be modified to run on the Center's computer. 
THe program was modified, and was Subsequently run using a check data 
set provided by Saxton. 


In addition to hydrologic factors which drive the model, such as 
daily precipitation, runoff or curve number, and pan evaporation, 
Saxton's model requires a knowledge of several crop or vegetation 
parameters or relationships, and soil characteristics. These 
include: the root distribution within the soil profile; the canopy 
cover with time during the year, and the Susceptibility of this cover 
to moisture stress; a vegetation phenology curve for the year, and 
the susceptibility of this curve to moisture stress 3) alicurve to 
represent the relationship between moisture availability and crop 
stress; the type, number of soil layers, and their thickness; and the 
dates when the range vegetation begins to grow and ceases to grow. 
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Some of the parameters, such as soil characteristics, may not be 
available, but can be fairly easily obtained by taking soil samples. 
However, vegetation relationships are not readily available in most 
range Situations. Therefore, much estimating is required, because 
the relationships presented by Saxton for agricultural crops in the 
Mid-west may not apply to rangeland vegetation on the Reynolds Creek 


Watershed. A sensitivity analysis, showing the effect of variations - 


in these relationships on predicted soil water, has not been 


performed to date, but will be important in determining what level of 
accuracy is required. 


After calibration of Wight's model using the 1979 Flats site data 
(see Interim Report No. 11), a continuous run for the years 1975 
through 1980 was performed. The measured soil moisture at the 
beginning of 1975 was used as a starting point, and then the 
predicted soil moisture at the end of each year was used as the 
beginning point for the next year. Results for each year are 
presented in Table 2.1. Shown are the predicted and observed soil 
moisture by layer, and the total water in the soil profile at the end 
of each year. The total water, as predicted by the model, differs 
from the observed amounts by up to nearly 30 percent. The greatest 
differences are associated with the wetter years, while predicted 
totals are almost identical to measured totals for the drier years of 
1979 and 1980. These results may indicate that the model as used 
removes too much water, thus reaching the lower limits for this soil 
nearly every year. Observation of predicted soil water in individual 
soil layers strengthens this suggestion. As shown in Table 2.1, the 
lower three layers reach essentially the same value at the end of 
each year, indicating that this value (which is different for each 
layer) is the lower limit of available water. 


Wight's model functions as a series of cascading reservoirs: thatsis,, 
all of the water infiltrated is assumed to be stored in the first 
layer until it reaches its field capacity, at which time any 
additional infiltration is assumed to be stored in the second layer, 
and so on. Since most of the roots are in the upper layers, and 
evaporation from the bare soil surface only involves the first layer, 
the majority of activity or water transfer takes place in the upper 
layers. Therefore, except in relatively wet years, the model would 
not show any water reaching the lower layers. This appears to be the 
situation shown in Table 2.1, where the first four years are 
relatively wet and measured water in the lower layers exceeds 
predicted. Whereas, for the last two drier years, the measured water 
and predicted water are essentially the same in the lower layers. 


An adequate comparision of the two models should include a comparison 
of factors needed to run the models in addition to final results. As 
noted previously (see Interim Report No. 11), these two models 
represent different levels of complexity. Saxton's model is more 
complex and attempts to include a representation of all physical 
processes involved, while Wight's model is a more simplified 
approach, which considers only broader or more general relationships. 
The data and model parameters needed also vary as shown in Table 2.2. 
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LAYER 
0 - 23cm 
23 - 46cm 
46 - 76cm 
76 - 106cm 
106 - 136cm 

TOTAL 
DIFFERENCE 


% DIFFERENCE 
















1975 1976 2 1978 1979 
WIGHT WIGHT WIGHT WIGHT WIGHT SAXTON 
PREDICTED MEASURED PREDICTED MEASURED PREDICTED MEASURED PREDICTED MEASURED PREDICTED MEASURED PREDICTED 
6.21 4.32 Bod 2.90 4.19 4.04 
















3.28 356 2.13 3.58 














2.68 





4.32 2.74 4.34 








3.56 






2.49 







4.83 





17.94 20.59 
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Table 2.1.--Comparison of Measured and Predicted soil water by layers using Wight's Model for 1975-1980 and Saxton's Model for 1979, in centimeters. 


1980 


WIGHT 
PREDICTED MEASURED 


BI ey) 4.06 
elt 3.18 
2.62 ao 3S 
2639 2.62 
32138 3.40 
16.23 16.59 
0.36cm 
2.2 


Table 2.2.--Input requirements for running Saxton's and Wight's soil water 


balance models. 


SAXTON'S MODEL 


inital sotl molsreurne 


Daily precipitation 
Daily runoff or curve number 


Daily pan evaporation 
Monthly pan coefficient 
Annual pan evaporation 


Soil type by layer 

Layer thickness 

Atime interval 

A soil pressure tolerance 


Root distribution 

Canopy cover curve 
Canopy susceptibility 
Phenology curve 
Phenology susceptibility 
Moisture-stress curves 
Planting & harvest dates 


WIGHT'S MODEL 


Initial soil moisture 


Daily precipitation 


Daily max & min temperature 
Daily solar radiation 


Soil moisture characteristics by layer 
Layer thickness 


Soil temperature curves by layer 
Canopy cover factor 
Plant growth curve 


Starting growth date 


Although most of the parameters needed can be obtained from the 
literature or data sets, with some modifications made for different 
site conditions, several of the parameters require information that 
is not readily available for rangeland sites. This is especially 
true for Saxton's model where six parameters or relationships 
concerning vegetation roots, growth, and stress are needed. In order 
to run the model, some information from the literature was used for 
root distribution, realizing that such data was collected from sites 
with different climatic and soil characteristics. Canopy and 
phenology curves were developed from vegetation measurements made at 
Reynolds Creek, also knowing that these data do not represent 


potential unstressed conditions. 


The canopy and phenology 


Susceptibility curves, and the moisture stress curves were used as 
presented by Saxton, even though these were developed from studies 
based on irrigated crops in the Mid-west. 


In using Wight's model, three factors require data that are not 
readily available. Of these, the plant growth curve, which is 
Similar to the phenology curve in Saxton's model, is the most 
difficult to make site specific, for the same reasons previously 
discussed. Data to develop soil temperature curves by layer are also 
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lacking; however, some models based on air temperature are available, 
and since the temperature of the deeper layers varies only slightly 
through the year, these curves can be estimated with some confidence. 
The third factor not readily available is the canopy cover factor, 
Which is merely a ratio of transpiring canopy to total canopy cover. 
Changes in this ratio of up to 50 percent made only minor changes in 
total soil water at years end; therefore, at least in this case, 
great accuracy was apparently not required in this factor. 


A comparison of the predicted and measured soil moisture, by layers, 
is presented in Figure 2.1 for both models. The ability of the 
models to predict variations in soil water status is of considerably 
more interest than the number of input items required, as long as 
these items can be obtained in reasonable time, and with reasonable 
accuracy. As noted in Figure 2.1, most of the changes in predicted 
and measured soil mositure occur in the top 46 em, or top two layers. 
In general, it appears that Wight's model follows the measured points 
somewhat better than Saxton's, except for the middle layer. The 
total amounts of moisture in the profile predicted by the two models 
are essentialy the same at the end of 1979, and bracket the measured 
value. Based on these results, it would be easy to conclude that 
Since both models produced essentially the same results, but Wight's 
model is easier to use, then why bother with the extra effort, time, 
and computer space required to use Saxton's model? However, these 
results are for only one year at a single location, and considerable 
effort was expended to make the predictions match the measured water 
content, since this was the calibration year. As noted in Table Cal 
Wight's model produced results that differed by nearly 30 percent 
from measured values for other years. Therefore, the tests will be 
continued on this and other sites to provide a more realistic 
comparison before conclusions are presented. 
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Figure 2.1.--Soil water content by layers predicted using Wight's model 
(solid line) and Saxton's model (dashed line) for 1979 at the 
Flats site, compared to measured values (represented by 
asterisks). 


It should also be noted that both models are still being modified, 
and the versions used in these tests are not the latest versions, but 
rather the ones in existence when the study began. In most cases, 
especially for Saxton's model, changes have been minor, and would 
probably not alter results significantly. The two major changes in 
Wight's model are the addition of a curve number based runoff 
routine, and the addition of built-in plant growth curves of various 
shapes that can be selected by the user to represent vegetation 
growth at a particular site. Since significant runoff events are 
rare at the Flats site, the runoff routine would not be likely to 
alter results. The effects of using the built-in growth curves 
ratner than the curve estimated from vegetation Surveys, are not 
known at this time. 
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Boise Front: A decision was made in April, 1981, between BLM and 
ARS, Northwest Watershed Research Center, to discontinue the 
vegetation surveys on the Boise Front rest-rotation System. This 
decision was made because an additional two years of data would not 
be enough to show the effects of the rest-rotation system, and State 
Fish and Game funding cuts precluded a project management level 
Suited for research purposes. It was agreed that the following 
research would substitute for the Boise Front study: 


Use Reynolds Creek vegetation and related resource data to 


evaluate rangeland monitoring techniques, key trend 
indicators, and sampling variability. 


Overstory cover sampling methods: Overstory cover was measured by 
the wheel-point (Tidmarsh and Havenga, 1915; and von Broembsen, 
1965), step-point (Evans and Love, 1957; and U.S. Department of 
Interior, 1979), and vertical point-frame (Goodall, 1952; National 
Academy of Sciences, 1962: and Hutchings and Pase, 1963) methods. 
Canopy cover was obtained from ten 100-point transects that were 
located in grazed areas at Flats (site 1), Nancy (site 2), and 
Whiskey Hill (site 3) study sites (see Figure 1). The wheel points 
were 2 feet from point to point. A pin guided by a notch in the boot 
was used in the step-point method. Plant hits with the wheel—point 
and step-point methods were recorded when any part of the pin hit a 
plant part. The point-frame was a 10-point frame with the pins 2.5 
inches apart and sharpened to a fine point on the tip: “Only hits on 
the point of the pin were recorded for the point-—frame method. 

For the three methods, hits were recorded from the ground up with a 
possibility of four hits per pin. For example, the first recorded 
hit (ground level) would be litter, rock, basal plant cover, or bare 
ground. Only rock .08 inches or more in diameter was counted. Moving 
up the pin, subsequent hits were recorded. Only the first three 
above-ground hits were recorded, and if the same Species was hit more 
than once per pin, only the first above-ground nit was recorded. In 
practice, the wheel-point and step-point pins were in place before 
the hits were observed; for the point-frame, the hits were observed 


as the pin was moved through the frame to ground level, but recorded 
from the ground up. 


Results from this study are presented in Tables een pand watt. =the 
results show that the wheel-—point and step-point methods give 
essentially the same results for both cover and frequency of 
occurrence. These results indicate that the wheel-point information 
obtained at the vegetation study sites on Reynolds Creek are directly 
related to survey data obtained by the step—point method used by BLM 
Survey crews. The point-frame method measured less cover and a lower 
frequency of occurrence of the grasses. These results would be 
expected because only tip hits were recorded by the point-—frame 
methcd and hits on the side of the pins were recorded for the other 
two methods. The pins were about 0.31 inches in diameter and 
represented a significantly larger area than did the pinsvips. 
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Table 2.3.-~Frequency of Bo mireaitece by the wheel-point, step-point, and point-frame methods, 














SITE 
1 2 3 

Wheel- Step- Point- Wheel- Step- Point- Wheel- Step-  Point- 

Point Point Frame Point Point Frame Point Point Frame 
Grasses 73 72 40 44 42 15 78 84 40 
Cheatgrass brome 65 66 35 we! i 0 63 70 33 
Sandberg bluegrass T T T 40 38 14 8 7 3 
Bottlebrush squirreltail 6 5 3 2 3 1 2 2 at 
Forbs 15 14 v 25 21 19 25 26 16 
Shrubs ¢ 14 14 8 18 16 17 41 35 27, 
Big sagebrush 7 7 4 18 15 16 27 23 18 


Shadscale 6 6 3 a a =a s— ae = 


l/ Frequency is the total number of hits divided by 1000 where no more than one hit per species was recorded for 
each point. 


2/ Trace. 
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Table 2.4.-Cover in Bon ae (includes foliar cover, mulch, rock, and bare ground) by the wheel-point, step-point, 
and point-frame methods. 














SITE 
1 2 3 
Whéeel- Step- Point- Wheel- Step- Point- Wheel- Step- Point- 
Point Point Frame ~~ Point Point Frame Point Point Frame 
Grasses 60 60 37 35 35 13 46 53 31 
Cheatgrass brome 56 57 33 T 1 0 35 44 26 
Sandberg bluegrass 72/ T T 33 : 33 12 6 6 2 
Bottlebrush squirreltail 3 2 3 1 1 1 1 1 1 
Forbs 7 7 5 10 8 13 8 7 9 
Shrubs 11 10 7 16 13 16 36 29 25 
Big sagebrush 6 6 4 15 13 16 25 2 18 
Shadscale 4 3 3 NA NA NA . NA NA NA 
Litter 3 2 16 6 5 13 4 4 25 
Rock 1 2 5 2 4 8 1 i i 
Bare ground . 18 19 30 31 35 37 rs 6 9 


i/ Cover is based on the first hit observed (going from the outer canopy to the ground surface) at each point, 
and is a percent of 1000 points for each method at each site. 
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Overall, the point-frame method measured grass cover as being 44 
percent less than the other methods. 


The hits of forb species were so few that the results were mixed, as 
was also the case for the shrubs at sites 1 and 2. At site 3 where 
the shrubs made up about 32 percent of the cover hits, the 
point-—frame method measured less than the other two methods by 23 
percent. As would be expected, the point-frame method measured more 
bare ground than the other methods. Because of the very small area 
of the pointed pin tip, it is less likely to encounter litter, rock, 
or basal plant cover. The higher cover of litter and rock measures 
by the point-frame method is somewhat unexpected, but may be due to 
the fact that as the wheel—point and step-—point pins were pushed into 
the ground, the litter and rock were pushed aside and weren't 
touching the pin when the observations were made. 


Point data obtained from the Flats and Nettleton sites were used to 
indicate how many step-points are required to measure basal cover 
within 10 and 20 percent of the population mean at the 80 and 90 
percent confidence levels. As expected, the data in Table 2.5 shows 
that the number of step-points varies according to the variable being 
sampled. If the variable sampled covers a high percent of the 
surface, fewer samples are required than when the variable covers 


only a small percent of the area, i.e.--the shrubs and forbs at the 
Flats site. 


Figure 2.2 shows the number of points required to measure cover at a 
given level of accuracy. As an example, if one wants to sample 
litter within 10 percent of the population mean at the 80 percent 
level of confidence, and litter covers 30 percent of the surface, 
line N2 shows that 400 points are required. 


The information in Table 2.6 was developed from basal cover data 
obtained from the Nettleton grazed and ungrazed sites in 1975. These 
data show that significant differences in cheatgrass basal cover 
could have been measured with relatively few points, because of the 
large difference in basal cover between the two areas. The data also 
shows that several thousand points would have been required to 
establish that the difference measured in the basal cover of shrubs 
was Statistically significant. In general, about 200 points in each 


area would be needed to detect cover differences of 5 percent or 
greater. 


Vegetation yield measurements: One of the procedures used by the BLM 
to measure vegetation yield is the weight-estimate method. This 
procedure was also used to measure the vegetation yield on the 
Reynolds Creek Watershed study areas. Data for two years from two 
Sites on Reynolds Creek Watershed were used to determine how many 
weight-estimates are required to estimate the mean yield, within a 
specific level of confidence (Table 2.7). 


The information in Table 2.7 indicates that 20 samples are required 
to measure all but three of the yield groups shown in the table to 
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Table 2.5.--Example of number of points required to measure basal cover, based on 1975 data 
from the Flats and Nettleton sites. 


—_-———————————— EES 
ee 


FLATS NETTLETON 








COVER CLASS MEASUREDE/ 2/ MEASURED 
COVER N1— N2 N3 N4 COVER N1 N2 N3 N4 


a 


Grasses 


Cheatgrass 141 1657 997 415— 250 ~259 867 522 217 130 

Total 2154 1497 900 374 225 -396 415 250 104 62 
Forbs - 026 9851 5926 2554 1536 -026 9851 5926 2554 1536 
Shrubs. -014 18800 11305 4699 2826 -030 8804 5296 2201 1324 
Total Live Cover »194 1132 681 283 170 451 Be ee 33) 199 83 50 
Litter -040 6535 3931 1634 983 -396 415 250 104 62 
Rock 129 1843 1109 459 276 044 6028 3626 1487 895 
Bare Ground -637 155 93 39 23 -109 2222 1337 557 335 


—_—_—__——— eee 
ee 


os Expressed as the ratio of actual hits to the total number possible (700). 

Z/ Nl = Number of points renntred to sample within 10% of the population mean at the 90% confidence level. 
N2 = Number of points required to sample within 10% of the population mean at the 80% confidence level. 
N3 = Number of points required to sample within 20% of the population mean at the 90% confidence level. 


N4 = Number of points required to sample within 20% of the population mean at the 80% confidence level. 


Number of Points 


20000 


1IOO0O Number of points required to 
sample within: 


Nl = 10% 
5000 90% 


N2 = 10% 
80% 


N3 = 10% 
90% 


N4 = 10% 
loo0 80% 


lOO 


30 


O lO 20 30 


of population mean at 
confidence level 


of population mean at 
confidence level 


of population mean at 
confidence level 


of population mean at 
confidence level 
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Figure 2.2.--Number of points required to measure cover in vegetation 


types represented by the Reynolds Creek Watershed. 
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Table 2.6,.--Example of number of points required to determine that the difference 
in basal cover between two areas was statistically significant, 
based on 1975 measured data from the Nettleton site. 


ee ee 


: 1 2 
Cover Basal taper 4 Number of Points” 
Class No Grazing Grazed N1 N2 


——————— a 


Grasses 


Cheatgrass 2239 053 19 eZ 
Sandberg bluegrass £123 piat 9238 5686 
Total - 396 jllele 26 16 
Forbs 026 -017 1401 858 
Shrubs - 030 -026 9068 5441 
Total live cover 451 235 25 16 
Litter ~ 390 293 114 70 
Rock -044 - 088 yal 104 
Bare ground .- 109 386 13 8 


eee 
1 


The ratio of actual hits to total hits possible (700). 


a2 Nl = Number of points needed in each study area to sample at the 90% 
level of confidence. 


N2 = Number of points needed in each Study area to sample at the 80% 
level of confidence. 
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Table 2.7.--Example of the number of weight estimates required to determine the vegetation yield 
(pounds/acre) at a sampling site. 


FLATS, 1977 PLALS eGo 
MEASURED SAMPLE sizet/ MEASURED SAMPLE SIZE 

MEAN YIELD N1 N2 N3 N4 MEAN YIELD N1 N2 N3 N4 

Grasses 
Bottlebrush Squirreltail 85 112 28 68 17 297 142 3D 86 2 
Total 101 72 18 43 ae) 720 93 23 56 14 
Shrubs 162 (MSc ais} 727/92 68 1324 129 32 77 19 
TOTAL YIELD 263 137 34 82 21 2061 61 15 37 9 

UPPER SHEEP, 1975 UPPER SHEEP, 1974 
MEASURED SAMPLE STZEL! MEASURED SAMPLE SIZE 

MEAN YIELD N1 N2 N3 N4 MEAN YIELD N1 N2 N3 N4 

Grasses 
Needle-and-thread 21 269 67 161 40 169 134 33 81 20 
Total 61 56 14 34 8 169 134 33 81 20 
Shrubs 764 164 41 99 as 1188 106 26 64 16 
TOTAL YIELD 1267 60 15 36 9 2105 32 8 19 5 


a Nl= number of weight-estimates required to sample within 10% of the population mean at the 90% 
confidence level. 
N2= number of weipht-estimates required to sample within 10% of the population mean at the 10 & 80% 
confidence level, 
N3= number of weight-estimates required to sample within 10% of the population mean at the 20 & 90% 
confidence level. 
N4= number of weight-estimates required to sample within 10% of the population mean at the 20 & 80% 


confidence level, 








within 20 percent of the population mean, at the 80 percent level of 
confidence. If the sampling program requires a mean yield estimate 
Within 10 percent of the population mean at the 90 percent level of 
confidence, about 250 weight-estimates would be required, and for 
Some years, up to 450 weight-estimates would have to be obtained. In 
general, the least weight-estimate samples are required for total 


yield estimates, and the greatest number are required “for shrubs and 
individual species. 


These results indicate the high variability of range vegetation 
parameters. To effectively monitor trend, inventory samples should 
be adequately large enough to detect ecologically significant 
differences. Indicator species should be relatively abundant, as 
well as ecologically important. The problem of monitoring trend is 
further confounded by the variability of climate, and better 
techniques need to be applied to factor out climatic effects. 


The wheel-point and step-point method produced similar results and 
could probably be used interchangeably. Both methods appear better 
adapted to routine trend analyses than the point-—frame method. 
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RUNOFF 


Plans programs and procedures; 
designs and constructs facilities 
for runoff studies; performs 
analyses and summarizes results. 


Streamflow and infiltration 
modeling. 


Runoff analyses and modeling. 


Tests various components in runoff 
models most applicable to 


Groundwater. 
Streamflow data collection, quality 


control, processing, and analyses. 


Performs data compilation and 
assists in analysis. 


Performs data compilation and 


processing. 


Designs, constructs, and services 
electronic sensors and radio 
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Snowmelt runoff. 


Snowmelt runoff. 
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1981 WATER YEAR RUNOFF CONDITIONS 


Watershed streamflows in the 1981 water year ranged from 53-102 
percent of the average of record. Snow accumulation on March Lise: 
1981 was near the lowest of record for that date because the 
October-February precipitation was only about 50 percent of average 
at 6800 feet elevation. However, precipitation the last halftof 
March and in April and May was 10-30 percent above average, resulting 
in peak streamflow rates in April at several stations. 


MICROWATERSHEDS 
Flats: No runoff was recorded at this station in 1981. 


Nancy: The only significant event was from 0.32-inch of rain 
February 13-14 when the soil was frozen, which produced less than 
0.01-inch of runoff from 3.1 acres. 





SOURCE WATERSHEDS 


Lower Sheep: Runoff from this 33-acre watershed was only 0.01-inch 


from 10 inches of precipitation, Table 3.1. Peak streamflow was 
about 0.1 ft /see on February 14. 


Reynolds Mountain East: Runoff from this 100-acre watershed at 6800 
feet elevation was 9.94 inches, about 50 percent of the 19-year 
average. Precipitation was 65 percent of average, 26.38 inches, 
Table 3.1. The stream was dry during September and most of August, 
Similar to the severe 1977 drought year. 


Reynolds Mountain West: Runoff from this 126-acre watershed was {gate 
inches, about 47 percent of average, the third lowest of record, 


Table 3.1. Peak streamflow from rain and snowmelt April 20 was about 
two-thirds of average. 


TRIBUTARY WATERSHEDS 


Salmon Creek: Runoff from this 8,900-acre watershed was 2.22 inches 
from 18.38 inches of precipitation, about 75 percent of average 
runoff, Table 3.2. Peak streamflow from nearly 2 inches of rain 
April 18-20 was about 40 percent of average. April runoff was about 
double the average of record. 


Macks Creek: Runoff from this 7,846-acre watershed was 1.46 inches, 


about 64 percent of average, Table 3.2. Peak streamflow was only 26 
percent of average, and occurred February 14 from 0.57-inch of 
precipitation associated with snowmelt. Peak streamflow from rain 
April 20 was nearly as great as in February. 


Dobson Creek: The station was discontinued in 1981. 
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Table 3.1. Water year precipitation, runoff, and peak streamflow, source watersheds, Reynolds Creek Watershed. 











z Reynolds Mountain East Watershed Reynolds Mountain West Watershed 


____Lower Sheep Watershed 











St 


Water Precipi- Peak Date Precipi- Peak Date Precipi- Peak Date of 
fe ee ec alo Runoff Streamflow of Peak tation Runoff Streamflow of Peak tation Runoff Streamflow of Peak 
inches inches £t>/sec inches inches ft/sec inches inches ft>/sec 

1963 16.98 a — — 37.82 11.11 4.16 Mprante eae ees ee 2 

1964 Soo <== --- --- 40.89 21.02 3.60 May 16 --- --- --- =- 

1965 20.86 --- --- --- 66.10 34.87 10.70 Dec. 23 --- 25.002/ 9.29 Dec. 23 
1966 6.81 --- --- --- 28.36 9.86 - 1.43 May 5 --- 7239 A cnehi/ Apia 13 
1967 18.73 0.342/ 1.41 Jan. 21 50.45 21.01 5.44 May 22 -- ial S 5.10 May 22 
1968 11.30 0.02 0.08 Feb. 18 ae OH, 6.72 1.48 Aug. 10 --- 6.31 Sey Feb. 23 
1969 1 32 0.52 0.49 Jan. 20 37.45 22.43 3.88 May 12 37 537 17.26 4.20 May 10 
1970 14.24 0.02 0.05 Jan. 27 39.60 20.06 5.89 May 17 37.95 20.24 12533 May 17 
1971 17.68 Oi Sl 0.20 Mar. 12 57.96 31.06 Sere May 4 45.75 21.41 10.24 May 4 
1972 13.82 0.91 2.08 Jan. 22 50.51 BSiney2 6.26 June 6 45.98 29.56 6.31 May 14 
1973 12.20 0.01 0.02 Apr. 17 31.01 13.24 Bost May 8 28.40 10.02 Sie o Apr. 27 
1974 10.28 0.26 0.38 Mar. 15 45.54 26.64 4.33 May 7 38.67 19.77 5.61 May 7 
1975 14.89 0.73 0.90 Feb. 13 Ly aes 27.93 O27 June 2 42.83 21.24 14,28 June 2 
1976 14.46 Ons 0.31 Mar. 17 42.51 22539 4.59 May 13 --- 16.38 4.09 May 2 
1977 8.27 0 0 --- Hu Alal 3.44 0.93 Apr. 16 == Pgs 0.72 Apr. 16 
1978 Tesi 0.14 0.09 Apr. 27 43.82 2302 4.50 May 14 --- 17.07 Se May 14 
1979 9.90 0.34 1.33 Jan. 11 32.42 1s) 53 BS} Bae May 15 --- 11.65 3.99 May 4 
1980 15.00 0.01 0.09 Jan. 12 41,37 18.79 3.69 Apr. 22 =—— 13.33 Any 22 Apia 22 
1981 10.02 0.01 0.10 Feb. 14 26.38 9.94 4.54 Apr. 20 - --- 7.28 3.83 Apree20 
MEAN E56) 0.30 0.53 = 40.89 19.59 4.59 ==: --- 15.49 5.70 ae 


First complete year of record. 


2 
=" precipitation record began in 1968 and terminated in 1975. 





Table 3.2. Water year precipitation runoff, and peak streamflow, Tributary Water- 
sheds, Reynolds Creek Experimental Watershed. 


_—_———O OO eee 


Salmon Creek Macks Creek 
Water Precipi- Peak Precipi— Peak 
Year tation Runoff Streamflow tation Runoff Streamflow 
inches | inches ft/sec inches inches ft?/sec 

1963 22.63 —— --- --- --- --- 
1964 SECS fo) --- --- --- --- =~ 
1965 S375 1 9.65 1007 --- =e 1200 
1966 10527 1.05 10 --- 0.61 12 
1967 PDL TED 2.24 85 ——< 1.54 90 
1968 14.73 0.77 34 --- 0.49 44 
1969 19.36 3.14 209 19.90 2298 307 
1970 24.96 3.07 210 19725 192 241 
1971 24.34 3.61 32, 23.65 Shs if, 281 
1972 22.74 5.50 201 PAD AB 4.84 138 
1973 OSES 2.14 55 UBISE NS 1.76 54 
1974 16.80 3.31 53 15.54 Bee 71 
1975 20.43 3.54 92 22.68 4.79 142 
1976 22.581 2.38 19 21.02 2.67 33 
1977 12.83 0.62 103 14.67 0.43 os TS 
1978 23.42 3.41 102 24.61 3.01 86 
1979 17.63 2529 380 16.45 231.0 300 
1980 Parana) Psa 43 20573 1.64 37 
1981 18. 38 Papier 69 15.69 1.46 47 


MEAS 20. 33 Sjsopl 165 Ue aenl 735 04i 182 


| 
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MAIN STEM WATERSHEDS 


Reynolds Creek Outlet: Runoff from this 57,700-acre watershed was 
1.40 inches, about 50 percent of the 19-year average, Table 3.3. 
Peak streamflow was about one-third of average from widespread rain 
April 18-20. Peak streamflow from rainfall and snowmelt February 
13-16 was slightly less than in April. The irrigation water supply 
was less than one-half of average in May and June, Table 3.4. 


Reynolds Creek Tollgate: Runoff from this 13,453-acre watershed was 
4.58 inches, about 57 percent of the 16-year mean, Table 3.3. Peak 
streamflow was 169 ft°/sec on February 16 from snowmelt and rain with 
frozen soil. April 20 peak streamflow was 149 ft°/sec. 

Precipitation was about 75 percent of the 19-year mean. Monthly 
runoff during May-July was only one-third of the average during this 
critical irrigation period, Table 3.4. However, crop yields in 
Reynolds valley were near or above normal because of favorable 
March-June rain. 


Table 3.3. Water year precipitation, runoff, and peak streamflow for main stem watersheds. 


Reynolds Creek Outlet 


Reynolds Creek at Tollgate 


Water ee Peak Date Precipi-~/ Peak Date 
tation Runoff Streamflow of Peak tation Runoff Streamflow of Peak 
inches inches ft?/sec inches inches ft/sec 

1963 25.03 1Batss) 2331 Janz 3 Sel -<~ --- --- 

1964 IBS 745, 2245 188 Jaleo 24.25 --- --- --- 

1965 26.83 Uses 3850 Dec. 23 38.93 ——= 11003/ --- 

1966 9.05 0.76 ao Apr. 1 ES TAS: $36.55) 5g Apr. 1 

1967 19.68 2.19 265 June 7 28.10 9.09 288 June 7 

1968 14.20 0.61 327 Feb. 21 pais pt 3.08 186 Feb. 21 

9 6 5aee 16.85 3.60 900 atime 29 alll 11.47 405 J ategys 

1970 20713 Pe oThe, 729 Jan. 27 ihe Bhs) 9.64 240 Sets, 227/ 

LOR 24.96 4.78 540 Jan. 18 41.89 14.98 196 May 

1972 eles 6.07 678 Mare 2 38 a2 16.45 Zi Mar. 

1973 16219 1.85 166 Apr. 17 25eL Be 6.00 147 AD emul /, 

1974 17.14 4.37 Zo Mar. 29 WSS WATS} 195 Mar. 29 

1975 LON, 4.12 281 Mar. 25 SUES is}. Sie 2a June 2 

1976 20.34 Pa Rath 140 Apr. 5 PES] S)G) 10.05 130 May 10 

ES, 11.42 0235 1119 June 11 15.49 ibeeyl 17 Apr. 8 

1978 19.64 Shs eas 589 Apr. 26 28.98 a Sz: 230 Apr. 26 

nono 14.56 2.06 1662 slequuer alah 20.19 6.78 121 Jans LL 

1980 18.51 Pan A Zoo May 6 29.08 S25 163 Apr. 23 

1981 14.40 1.40 251 Apr. 20 AAD eAS)ah 4.58 169 Feb. 16 

MEAN Veaie2 5 2.87 770 Iba] CN) 8.91 191 


Raingage No. 116X91. 


Ay 
2/ Rainzgage No. 155X07. 
3/ Estimated peak flow. 


| 





Table 3.4. Water year runoff in 1981 and the mean of record by months. 











Reynolds Creek Reynolds Creek 
Month Outlet Runoff Tollgate Runoff 

IO) SAO 31 on. ole 0b LoG lL 

ME oa Fam as Se lames GY ORV cae inches ~———-————+ => 
October AAS) 02025 076 0.082 
November - 036 0.046 083 Do 2) 
December «99 Oe 61 2 2S 
January | OVS 0.380 a WAS) eh 
February poi) Oe 2 0.422 
March 145 0.455 ~405 O 7 265 
April Pas 0.550 16/3 1.788 
May 248 0.609 a Wha lls 3.088 
June «20 Oee2.96 ~3%3 12358 
July 024 0.048 050 0, 234 
August G12 Oe0 27 2OUu 0.046 
September 003 02013 _-001 0.037 
Total 1.403 22070 4.581 Bh9L0 
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INFILTRATION BASED RUNOFF ESTIMATION 


INTRODUCTION 


When compared with the SCS curve number procedure for runoff 
prediction, the infiltration approach requires significantly more 
input data and a much more complex computational methodology. 
However, the general availability of soil and soil water data is 
increasing the feasibility of the infiltration approach. The 
infiltration approach also enhances our ability to evaluate the 
effects that land use practices and treatments have on the generation 
Of PunoLs s 


Table 3.5 compares the factors in the SCS runoff equation procedure 
with the infiltration procedure. Factors correspond to those that 
are described in the SCS National Engineering Handbook, Section 
4-Hydrology, Tables 4.2 and 9.1, and Fig. 10 (1972). The 
Precipitation Factor required in the infiltration approach is a 
rainfall distribution and amount, as compared to the simple amount. 
The Soil Factor under the infiltration approach represents the 
initial soil-water condition and the soil infiltration equation 
parameters. Cover Factors represent the influence that land use and 
soil surface practices have on the infiltration equation. This 
influence is assumed to relate to soil hydraulic properties, such as 
Soil porosity, and the antecedent water content. Surface storages, 
surface interception and surface depressions, are determined by 


surface cover and surface practices, as discussed by Dunne and 
Leopold (1978). 


The infiltration approach presented here applies the Green and Ampt 
infiltration equation to the SCS Type II, 24-hour rainfall 
distribution. Table 3.6 presents the 24-hour rainfall intensity and 
accumulation distributions, respectively (Kent, 1968). 


Table 3.7 presents the estimates for the Green and Ampt parameters. 
These values were developed in studies by Rawls et al. (1981). The 
infiltration approach must be applied in two Stages during the 
rainfall event, i.e., pre-ponding and post-ponding. Prior to surface 
ponding, the infiltration rate is the rainfall rate minus the surface 
interception rate and, thus, is generally rainfall dependent. After 
surface ponding, the infiltration rate is determined by soit 
properties as long as rainfall excess continues. The two Stages are 


linked by a time correction, which is applied to the rainfall time 
scale. 


oe 





Table 3.5.--Comparison of direct runoff prediction approaches. 





Factors 


Approach 





1/ 


Curve Number— 


Infiltration 


_--—_———————————————————————————— ee 


Precipitation: 


Soil 


Cover 


Storage: 


Soil 


Surface 


Interception 


Rainfall amount 
(24-hr) 


Antecedent Moisture 
Gondition) (isn hi or Tin) 
(5-day prior rainfall) 


Hydrologic Soil Group 
(HSG) 
CAS Bis GS com. D) 


Land use 
Treatment or Practice 
Hydrologic condition 


Initial abstraction (I,), 
assumed as 0.2 (S) 


Included with initial 


abstraction, Ig 


Included with initial 
abstraction, I, 


Rainfall intensities or 
rainfall distribution and 
amount 


Antecedent soil-water storage 
(volume) by soil layers or 
soil depth 


Soil-water properties by 
layers or soil depth, i.e., 
bulk density, saturated 
conductivity, and soil- 
water entry or bubbling 
pressure 


Influences on soil properties 
by: 

Tillage 

Amendments 

Traffic/animal compaction 


Crusts, desert pavement 
Mulch cover 


Vegetation cover 


Infiltration prior to surface 
ponding 


Estimated surface depressional 
Storage as influenced by 
topography and land use 


Estimated interception storage 
by ground cover (live and/or 
mulch) 


__ SSSSSSSSSSSSSSSSSSSSSSsSSSssSsSsSsSsssesssSsSsseF 


1/ 


— The factors under the Curve Number Approach are taken from Table 9.1 and 
Figure 10.1, SCS National Engineering Handbook, Section 4 - Hydrology. 
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Table 3.6.--SCS Type II, 24-hour duration rainfall distribution. 





Time Accumulated Rainfall Time Accumulated Raincal = 
rainfall amount/ intensity/ rainfall amount/ intensity/ 
24-hr total 24-hr total 24-hr total 24-nr total 








10.0 
1025 
11.0 
ONS S. 


4l 








Table 3.7.--Green and Ampt infiltration equation parameters - mean values 
for each soil texture class. 








Soil Total Effective Wetted Front Conductivity 
Texture Porosity Porosity Capillary 
Pressure 
PO WF K 

(=) (>) (cm) (cm/hr) 
Sand 0.438 0.417 4.95 11.78 
Loamy 
Sand -437 -401 6.13 2.99 
Sandy 
Loam -453 ~412 11.01 1.09 
Loam -463 434 8.89 - 340 
Silt 
Loam 501 486 16.68 ; -648 
Sandy 2 
Clay Loam - 398 - 330 21.85 +153 
Clay 5 
Loam -464 390 20.88 -097 
Silty 
Clay Loam -471 432 27.30 097 
Sandy 
Clay -430 321 23.90 - 064 
Silty 
Clay 479 423 Po) ae -O51 
Clay 475 - 385 31.63 034 
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In the following sections, the procedures are developed and, by 
examples, are illustrated for predicting the surface runoff amount 
for the 24-hr, SCS Type II rainfall distribution with the Green and 
Ampt infiltration equation. The specific computational procedures 
are available in a BASIC computer program written for an APPLE II 
Plus with 48 K RAM and one DIS es The basis of the procedure is an 
iterative solution of the developed equations with the tabulated 
rainfall distribution values in Table 3.6. Estimates of peak runoff 
rates can be made from the predicted runoff volume and duration 
utilizing the procedure in Kent (1968), or the equation developed by 
Smith and Williams (1980). 

PRE-PONDING INFILTRATION 
Time to ponding 
At the time of surface ponding, the infiltration rate equals the 
rainfall rate. The amount of infiltration up to that time is the 


accumulated rainfall minus interception storage. At that time the 
Green and Ampt infiltration equation is written as: 


f = p = K[1 + nWF/(P-SI)] [1] 


where 
f = infiltration rate, em/hr 
p = rainfall intensity, ecm/hr 
P = accumulated rainfall, em 
SI = surface interception, ecm 
K = conductivity, em/hr (Table 3.7) 


nN = available soil porosity, effective porosity reduced by 
antecedent soil-water = (PO-AW) 


WF = matric pressure at the wetting front, cm (Table 3.7) 
PO = effective porosity (Table 3.7) 


AW = antecedent soil water 


Trade names and company names are included for the benefit of the 
reader and do not imply any endorsement or preferential treatment of 
the product listed by the U.S. Department of Agriculture. 
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Since the rainfall rates and amounts are tabulated (Table 3.6) as a 


ratio to the 24-hour amount, P2, equation [1] is re-written as: 


OM 
P2CP Ly seeker lL bw P2(PA-SI/P2) [2] 


where 
OM = WF(PO-AW) 


PI = p/P2, hr~!, the intensity ratio 


PA = P/P2 
P2t=.total <24—-hr rainfall?) cm 
OP = effective porosity 


SI = surface interception 


Multiplying both sides of [2] by P2/K(0M) gives: 


(P2) aP2 i, 1 
Keay)? = OM * (PRCSI7P2) 


or 


2 ee 
Ws) latPt)1— 8S + Caamenypapaderee [3] 


where 


PS P2/0M 


R 
ul 


OM/K 
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Solving for PS gives 


PS = 1+ V1+ 4 (PI)/(PA-SI/P2) 


oewoers case eae C4] 


Equation [4] is solved by an iterative procedure using Table 3.6. 
The time at which the Right Hand Side (RHS) of equation [4] equals 
the known Left Hand Side (LHS) of [4] is the time the infiltration 
rate equals the rainfall rate. Inputs are rainfall amounts and 
intensity ratios, PA and PI, the soil texture, the antecedent soil 
water content, and the 24-hour rainfall amount. Computation 
restrictions placed on SI and DS are: 


nism so6(P2) 


DSLcere 


Time correction 


After surface ponding, the infiltration rate proceeds at the maximum 
rate which the soil will sustain until the rainfall rate falls below 
the infiltration rate. The integrated form of the Green and Ampt 
equation then applies: 


FC — (OM)£n(1 + =) = KT [5] 


Where FC is the accumulated infiltration amount, T is time, and (OM) 
and K are the parameters defined for equation [1]. However, the time 
Scale in [5] must be adjusted for the pre-ponding time period, i.e., 
initialized to the start of rainfall. The time correction was 


developed by Mein and Larson (1971). The correction for the time 
Seale is: 


TC = @OT BA OTE) [6] 


where TP is the time to ponding and TE is the equivalent time for the 
pre-ponding infiltration if it had occurred with ponded conditions. 
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Calculation of TE is facilitated by an approximation of equation [5] 
derived by Li et al. (1976): 


TE. = (F1)2/1K(2(0M). + F1)] {7] 


Where F1 is the infiltration amount at ponding exclusive of surface 
interception. 


POST-PONDING INFILTRATION 


Infiltration rates and amounts 


After ponding the approximation derived by, Lt et al... (1976). is-used 
to calculate accumulated infiltration 


Pio ON / eet Lobel LL O8 + TT. [8a] 


where, TT: = (T.~ TC)/a 


R 


OM/K 


The infiltration rate is calculated as: 
Poy Gl ot OM EC) [8b] 
where PF.=. infiltration; rate; cm/hr 
Cessation of Rainfall Excess 
At the termination of rainfall excess, the infiltration rate again 


enuats cone rainfall rate: 


p = K(1 + OM/FC) 
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or 


1 


P2 ; Pita 
K (Uy mae (FC/OM) 


or 


PS (aPI) = 1 + 1/FS [9] 


where FS = FC/OM 
PSi = P2/K 


PI = p/P2 


Solving [9] for FS gives: 


FOog= 17 Po, (oP Leb [10] 


Equation [10] represents the total amount of infiltration that 
reduces the infiltration rate to the appropriate rainfall rate. 
Substituting equation [10] into the Li approximation (equation [7]) 
and solving for PS gives, after simplification: 


bho bn / 8a 
PS = Wo (PI) 3 + ies TTC (alsa 


where a = OM/K 
Time in [11] is limited to the interval: 


T2c<e Ts. 2h 


Equation [11] is also used with an iterative procedure and Table 3.6 
to interpolate the time at which the rainfall rate equals the 
infiltration rate, i.e., when the right hand side equals the left 
hand side. 
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DISCUSSION 


The application of the runoff procedures is illustrated with the 
examples presented in the following section. In Figures 3.1 and Jae. 
the flow charts present principal computational Segments of the 
program. The SCS Type II rainfall distribution resides in the 
Program; however, other 24-hour distributions could be substituted. 
Soil inputs are given in Table 3.7 by soil texture. The 24-hour 
rainfall total is program input. 


Equation [4] is the basis for calculating the time of ponding, TP. 
In principle, values of accumulated rainfall ratios, PA, and rainfall 
intensity ratios, PI, for succession times, T, in Table 3.6 are 
Substituted in the Right Hand Side (RHS) of equation [4] and its 
values compared with the known Left Hand Side (LHS). The lower and 
upper bound times are established, which bracket the LHS of equation 
[4]. A linear interpolation equation is then used to calculate the 
ponding time. If surface ponding does not occur by the end of the 
initial 12 hours, then no rainfall excess occurs. Note that surface 
interception, SI, is excluded from rainfall was itAdoes not 
contribute to soil infiltration. From the ponding time, the 
infiltration at ponding can be calculated. The time correction, 
equation [6], is calculated. Soil surface storage filled during the 
initial 12 hours is also calculated. Figure 3.3 shows the time 
positioning of these calculated quantities for the first example. 
Calculations for a second example are presented to show a situation 
where depressional storage is not filled in the initial 12 hours. 


Following surface ponding and during the period of rainfall excess, 
infiltration amounts can be calculated using equation [8a] and rates 


by equation [8b]. The rate curve is plotted in Fig. 3.3 Staring at 
ponding. 


The second half of the program uses the final 12 hours of the 
rainfall distribution input (Table 3.6). Equation [11] is the basis 
for calculating the cessation of rainfall excess. By substituting 
Sucessive values of PI and Time from Table 3.6, a lower and upper 
bound is established whereby the Right Hand Side of [11] brackets the 
known Left Hand Side of [11]. Again, a linear interpolation equation 
calculates the Time that rainfall excess ends,.TX....The total 
rainfall excess is then calculated and soil surface Storage is filled 
by available rainfall excess. If soil surface storage is filled, 
then surface runoff is calculated and its duration is determined. If 
the surface storage component is not filled by rainfall excess, then 
no surface runoff is calculated. 


Tne output listing for the examples show the calculated quantities. 
Figure 3.3 shows the corresponding time sequence of the calculated 
quantities for example 1. Table 3.8 presents a list of symbols used 
for input, computations, and output. 
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Sott - K, PO, WF 
SOIL WATER - AW 
SToraces - SI, DS 
24-HR RAIN - P2 

















CALCULATE: 
DIFFERENCE BETWEEN 

RHS anp LHS oF Ea, 4 
AT SUCCESSIVE TIMES 










RAINFALL 
TABLE 
EXHAUSTED? 


New 
LOWER BOUND 
on TP? 






NES Qutput: 


No RUNOFF 






UPPER 
BOUND ON IP? 


YES 


CALCULATE: 
PONDING TIME 


INFILTRATION AT 
PONDING 


TIME CORRECTION 


CALCULATE: 
DS FILLED IN FIRST 
12 HouRS 





Pigurerd. lini tials hours ote tata iL. 
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Lyput: FINAL 
12-HouR RAINFALL 
DISTRIBUTION 










YIFFERENCE BETWEEN 
RHS and LS oF 

ed, il at 
SUCCESSIVE TIMES 






NEW 
UPPER BOUND 
on TX? 







RAINFALL 
TABLE 
EXHAUSTED? 










YES 














LowER BOUND 
on TX? 


TX = 24 yours 


CaccuLate TX 






CALCULATE RAINFALL 
Excess, RX 










SURFACE 
STORAGE 
FILLED? 





No RUNOFF 







CALCULATE RUNOFF, 
Q, RUNOFF DURATION, 
TQ 






———————eee 


| 


Figure 3.2.--Final 12 hours of rainfall. 
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(PA = P/P2) 


(Qi ah 
Og 
- 08 


-07 


-06 


-05 


-04 


-03 


-02 


sO k 


Figure 3./3.--Example 1. 


CLAY LOAM 

AW = 0.25 cm 
Sia aC acm 
DSe=, 3. 0aeem 
P2 = 20.0 cm 
FP = 0.96 cm 
TP = 6.98 hr 

Q = 11.75 cm 
TOR= si 2eoleanty 

Q 
12 V5 18 21 


TIME (hours) 


JL 


24 





—_— EEE 





EXAMPLES 


CHLaoOLEwTEXTURE }-6Clay~Lean® 


INPUTS: Ke =0. 097° cm/hr , 24 AWrS Ores, eS POMa sg. 
Hiw=r2Os Oo. CM eo. = 1. 0 ome aloe 3.0 em, 
P2 = 20 am 
QULBUTS: 
PPS “4¥964em TP = 6.98 hr 
Fil = -96 cm Diese -20Gbrirs 
Sis 12 02 ein See Bee Oe hrs 
DS = 3 ROM TD =° 17379Snrs 
RX =) 20.0 em TX = 2430" hrs FT = 4.246 em 
Runoff is produced: 
Qe 94 75 Jom Wes) eeo | ane 


C2)ee SOIL TEXTURE: Sand 


INPUTS: K = 11.688 ecm/hr, AW = 0 
Sve 1. 0%ems (Ds =" 5. 0-em 


OUTPUTS: 
FP = 16.93 cm EPPS S87 
Baie= 155 938 %en TASSe B98 
Slest 430 “om Slee oOo 
Ds <0 h5.0°-em TDF 22 at 
RX = 20.49 cm TX = 2.14 


Runoff is not produced: 


DS (FILLED) 


aes 


roo, UTS tly, We = 4.05 om, 
were 30S eM 


hr 
hr 


hrs 
nours 


hr FI*= 79.08 om 


a= 0541 cn 


Table 3.8.--Input and output parameters. 


Rainfall Distribution Input (Table 3.6) 


PA Tabular accumulated rainfall amounts divided by 24-hour total - 


PI Tabular rainfall intensities divided by 24-hour total, hr-1 


T Tabular time, hr 


Soil, Storage and Storm Input 


K Green-—Ampt conductivity, em/hr 
AW Antecedent water content, em? /om? 


PO Effective porosity, total porosity minus residual saturation, 


WF 


SI 


DS 


Fe 


Output 


FFP 


Fa 
ZR 
nl 


SI 
ST 


DS 
TD 


RX 


TX 
a 


TQ 





em?/em? 


Green—Ampt wetting front suction, em 


Surface interception, cm 
(must be less than 66% of 24-hour rainfall) 


Soil surface storage, cm 
(must be less than 24-hour rainfall) 


24—nour rainfall total, om 


Infiltration total at surface ponding, including surface 
interception, om 


BP minus: o2 
Time at surface ponding, hrs 
TP minus time to fill Si 


Surface interception, em 
Time at which surface interception is filled, hr 


Surface storage, em 
Time at which surface storage is filled, hr 


Total rainfall excess, cm 
Duration of RX, hr 

Total infiltration, em 
Surface runoff volume, cm 
Duratian of 0, hr 


2S, 











CONCLUSIONS 


Equations are presented for applying the Green and Ampt infiltration 
equation to an SCS Type II (24-hour) rainfall distribution. Other 
24-hour distributions, as well as a constant intensity distribution, 
can be placed in the program and runoff calculated. Soil inputs for 
the infiltration parameters are included. The predicted runoff 
amounts and durations can be utilized in making peak rate estimates. 


For design situations, the infiltration approach provides an improved 
framework for evaluating the effects of land use and practices on 
runoff amounts and duration. For example, changes in soil porosity 
can be directly entered as an infiltration equation parameter. The 
availability and ease of operation of micro-—computers enables the 


practicing hydrologist to utilize runoff models, such as the one 
reported here. 
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Reynolds Creek (Reynolds Creek Experimental Watershed station 
locations are shown in the Introduction, Figure 1). 


MICROWATERSHEDS 
Flats: No runoff or sediment yield in 1981. 


Nancy Gulch: Insignificant sediment yield in 1981. 


SOURCE WATERSHED 


Reynolds Mountain East: Total sediment yield from this 100-acre 


watershed in 1981 was 10.1 tons, 78 percent of the 14-year mean, 


Table 4.1. The maximum suspended sediment concentration was 565 mg/1 
with peak streamflow. 


TRIBUTARY WATERSHED 


Macks Creek: Suspended sediment yield from this 7,846-acre watershed 


was 143 tons, only 7 percent of the 14-year mean, Table 4.1. Nearly 


70 percent of 1981 sediment yield was on February 14, and April 20. 
Tne maximum suspended sediment concentration was 1,990 mg/l on 
February 14. 


MAIN STEM WATERSHEDS 


Reynolds Creek at Outlet: Suspended sediment yield from the 
57,754-acre watershed was 1,867 tons in 1981, 15 percent of the 
15-year mean, Table 4.1. About 88 percent of the sediment yield was 
from snowmelt associated runoff in February and April. The maximum 
Suspended sediment concentration was about 6,000 mg/l on April 20. 


Reynolds Creek at Tollgate: Total sediment yield from this 
13,453-acre watershed was 1,340 tons, 27 percent of the 15-year mean, 
Table 4.1. About 85 percent of the sediment yield was produced in 
peak runoff periods in February and April. The maximum suspended 


sediment concentration was about 6,000 mg/l with peak streamflow 
February 16. 


VEGETATION GROUND COVER ANALYSIS AND PREDICTION 


Percent bare ground at peak standing crop, as an indicator of erosion 
potential, ranged widely from year-to-year on vegetation study sites 
on the Reynolds Creek Experimental Watershed, Table 4.2. The 
objective of this study was to predict on April 1 of each year, the 
percent bare ground at peak standing crop based on: 
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Table 4.1.--Sediment yield in tons at Reynolds Creek Watershed Stations. 





ea DERE ES BEES SEE ESI PTI Py a ee eT ee am fe 











Reynolds 1/ Reynolds Reynolds— 

Mountain Macks— Creek at Creek at 
Year East Creek Tollgate Outlet 
1967 -——— --- VBS) 13503 
1968 es) 393 1965 4334 
1969 17.0 6332 12994 39336 
1970 splat 3585 7242 15369 
iRrede 18.1 5833 O77 28641 
SY 18.3 5414 8838 37396 
1973 9.4 1147 1203 2415 
1974 L053 1214 2774 5762 
LOTS 14.2 1949 7867 9860 
197.6 12.4 646 2546 1430 
1977 1.0 7 51 3257 
1978 225) 554 2797 8256 
1979 Shere 1634 1808 11674 
1980 Wars 87 1778 4237 
1981 LORE 143 1340 1867 
MEAN TS20 2067 4950 12489 








1 
as Suspended sediment only. 


Table 4.2.--Percent bare ground at peak standing crop (PSC), Reynolds Creek 
Watershed vegetation study sites. 





Site 1972 oS. 1974 1975 1976 METE bs SESIIKS: NES 1980 


eae rare percent bare ground---~---------------~--—--~------ 


057 Ungrazed 40.6 54.6 34.8 69.2 89.7 44.8 55.6 39.4 oe oot 


Grazed 39.3) 54.3 39.0 63.4 77.0 62.5 Whe st 6l.1 42.6 
092 Ungrazed -- 37.0 26.6 26.8 31.8 SIS ate! CHESS) 29.8 1224 
Grazed -- 29.7 33.0 29.4 23.9 34.6 37.6 37.0 24.2 
098 Ungrazed She Hf 54.6 32.9 38.2 23.6 30.5 46.2 40.5 22.0 
Grazed 21.9 54.3 says 48.5 eke eos 0 47.3 43.1 29.3 
135 Ungrazed 16.7 24.8 17.9 11.9 12.0 24500 wlaas 28 10.4 
Grazed 16.9 Sistats: 30.6 Bond 2159 33.0 33.0 39.1 Zine 
117 Ungrazed 16.8 -- 26.8 25.07 Ba 2 Oren 14.5 aoe Vie 
Grazed 15.4 -- Pipes: 3128 6.9 38.1 33.0 Boo L526 
138N Ungrazed 5.4 14.3 £155 6.9 V8 15.0 Ses 3.0 dee 
Grazed 9.53 18.0 Le 5 6 2220 es 18.0 17.4 PaO 
138S Ungrazed 18.1 -- BL AG) 40.1 36.4 Saleay/ 40.0 47.8 3254 
Grazed 26.9 -- 48.1 47.0 26.6 65.0 41.4 54.8 41.4 
176E Ungrazed Lig Lek 24.6 19.6 3H be(O) 1 ets 1536 21.8 18.8 
Grazed 13.9 19:58 Zou 23.6 29.4 21.8 16.6 2760 29.2 
176W Ungrazed Pe 2 -- LT oak 24.4 14.3 5.) agi 8 Lot Teo 
Grazed 19.9 -- 2aeT 29°.9 6.4 46.3 16.7 PBC) Dae 
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1. The percent bare ground at peak standing crop the previous 
year. 


2. The percent of the area not protected by overstory vegetation 
at peak standing crop the previous year. 


3. The accumulated precipitation in inches at the site during the 
previous 6-month period, October—March. 


4, Estimated precipitation in inches during the April-—July 
growing season of the current year. 


A precipitation analysis showed no significant persistence between 
precipitation amounts in the October-—March winter precipitation 
accumulation period, and amounts in the following April—July period 
when most vegetative growth occurs. Consequently, no estimates were 
made of April—July precipitation. 


The generalized form of the multiple regression equation used in the 
analysis of bare ground data was: 


Y = Do + by X, + boXo + D3X3 


where Y is the dependent variable, percent bare ground at peak 
standing crop of the current year; X, is an independent variable, 
percent bare ground at peak standing crop the previous year; X2 is an 
independent variable, the percent of the area without overstory 
canopy; and X; is an independent variable, the accumulated 
precipitation during the previous 6-month period October-—March; and 
Do» by, bo, and bg are regression coefficients. Regression equations 
developed in the analysis are summarized in Table 4.3. The 
coefficients of multiple regression are a measure of the degree of 
association between observed and predicted values of percent bare 
ground at each site, and results range from excellent to poor for the 
period of study. Overall, predictions of percent bare ground are 
better for ungrazed areas. 


Further improvement in April 1 predictions of percent bare ground, Y, 
was attempted uSing percent bare ground at the end of the grazing 
season, instead of at peak standing crop, X,;. Equations and 
correlation coefficients from this analysis are shown in Table 4.4. 
Since data were collected at the end of the grazing season for only 
four years of the 9-year study period, equations were not developed 
for individual sites. Overall, use of data collected at the end of 
the grazing season slightly improved the predictions. 


The procedure used in this study has application in managing 
rangelands through drought and wet years when potential erosion 
ranges widely. Also, the procedure could be useful in trend plot 
analysis for sites with a good historical record of bare ground data. 
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Table 4.3.--Regression equations developed for individual vegetation sites 
on the Reynolds Creek Experimental Watershed, 1972-80 data at 
peak standing crop (PSC). 


a 





1 

Site Equation ay 
057. Ungrazed Y = -36.30 + 0.70X, - 0.49X> + 13.59X3 0.66 
057 Grazed Y= 10.25 + 0.64%, - 0.34X5+ 5.15X3 0.46 
092 Ungrazed Y= 20.32 + 0548xX) =10.68X%5 —  0202K3 0.47 
092 Grazed Y=) 35. 16greOrI 2X ete 05 07 Xioa- nO a7 oxy 0.49 
098 Ungrazed Y= 40.14 - 0.62X; + 0.44X9 + 1.78X3 0.44 
098 Grazed Y = (53.19 — 2 0x 0.62% oe 2531k3 0.81 
117 Ungrazed Y = - 36.97 = 1.03X)/+.0.70Xo 4 0702X3 0.73 
117 Grazed ¥Y = 44.90 — 0.25xX) + 0.86%5 — 9 3.65X3 0.68 
135 Ungrazed NE 10.44 + 1.47X, - 0.82X2- 0.90K3 OF7u 
135 Grazed Y= 845.20 — 0738x, + 005k 5 — 0211 xX5 0.54 
138N Ungrazed Y= 15.18 - 0.31X; + 2.16X5- 1.04K, 0.73 
138N Grazed Views (LOSG49-+00.2 7h — el 09k oem Oo sca 0.70 
138S Ungrazed Y= 48.69 + 0.50X; - 0.84X5 = 1.19X3 0.97 
138S Grazed Y= 80.52 - 0.53X; + 0.84X5 - 3.68K3 0599 
176E Ungrazed Y =) 20.02"— 0747) tel s23x%> + O0.05X3 0.67 
176E Grazed Y= 15.48 + 0.24X, - 0.04K,+ 0.12X3 0533 
176W Ungrazed Y= 60.56 —"1. 66%) “+13 Sika 2. xX 0.70 
176W Grazed Ve= D4. SD e ae le C4 ky ere ON oe Uns 0.80 
ALL Ungrazed Y= 14.40 + 0.43X; + 0.34X9 - 0.36X3 0.70 
ALL Grazed Y= 24.88 + 0.34%, + O0.17X9 - 0.51X3 0.62 
SE Da ag 
1/ 


— Coefficient of multiple correlation. 


Table 4.4.--Regression equation developed from data collected at the end 
of the grazing season (EGS) compared with data collected at 
peak standing crop (PSC); data from 1974, 1975, 1978, and 
1979 used in both analyses. 








Data Base Equation ie 

PSC Ungrazed N26 Pk USO Kk ne sOOO Ns 0.76 
PSC Grazed Y = 3.34 + 0.68X,; + 0.20X5 + 0.04X3 0.69 
EGS Ungrazed = 11.34 + 1.64X, - 0.96% - 0.57X3 ? 0.78 
EGS Grazed =a e 2d tel OUR a=) Ooo Xie s O eee a 0.74 
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SOIL LOSS STUDIES 


The Universal Soil Loss, Equation (USLE) tis 
Ay = RKLSCP 


where A is the computed soil loss per unit area; R is the rainfall 
and grunort facton; eKeisptheysoil serodibility factor: .L<is ‘the 
Slope-length factor; S is the slope-steepness factor; C is the cover 
and management factor; and P is the support practice factor 
(Wischmeier and Smith, 1978). Application of this equation to 
prediction of soil losses from sagebrush rangelands is poorly 
understood, because basic data is lacking on USLE factor values for 
rangeland conditions. Problems associated with applying the USLE to 
the Reynolds Creek Watershed are: 


1. The rainfall and runoff factor, R, does not properly account 
for snowmelt and frozen soil runoff. 


2. Soil erodibility, K, values for rock pavements and rocky soils 
are questionable. 


3. Slope-length factor, L, values on sagebrush rangelands are 
probably much less than on cultivated lands, because rills are 
usually short and discontinuous on rangeland slopes. 


4, Cover and management factor, C, values for rangelands with 
Well distributed dense-rooted grasses and shrubs are commonly 
overestimated. Rangeland vegetation is probably more 
effective in controlling soil loss than evenly spaced 


cultivated crops with a similar percent ground and canopy 
cover. 


Presently, the greatest need in applying the USLE to rangelands is 
developing a procedure similar to that of Dissmeyer and Foster, 1980, 
for forest lands. To accomplish this, standard rainfall simulation 
erosion measurements and natural runoff plots are needed to provide 
basic data, similar to extensive data used in development of the USLE 
for cultivated areas. Meanwhile, the USLE is useful in predicting 
relative effects of grazing and brush control on potential erosion 
(Johnson, et al., 1980; Northwest Watershed Research Center, 1980). 


SEDIMENT YIELD STUDIES 


Pacific Southwest Inter-agency Committee (PSIAC) procedure: 

The PSIAC procedure (Pacific Southwest Inter-agency Committee, 1968) 
was adapted for using Reynolds Creek Watershed data (Northwest 
Watershed Research Center, 1979). Predicted sediment yields were 
Within about 15 percent of measured watershed sediment yields. The 
PSIAC ground cover, land use, and upland erosion factors increased 
with grazing and predicted sediment yields were about 14 percent 
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greater on grazed than on ungrazed areas. Sagebrush eradication by 
Spraying increased predicted sediment yields about 8 percent, and 
brush eradication by cutting and removal increased predicted yields 
about 12 percent on ungrazed areas (Johnson and Gebhardt, 1982). 


Results of this study show the PSIAC procedure can be used to predict 
long-term watershed sediment yields from sagebrusn rangelands where 
data are available to evaluate the factors. Also, the PSIAC factors 


influenced by changes in land use appear to respond properly for the 
eases tested. 


Modified Universal Soil Loss Equation (MUSLE): 


The MUSLE was developed by Williams, 1975, to eliminate the sediment 


delivery ratio when applying the USLE to watersheds and specific 
storm runoff and sediment data. Williams replaced the conventional 
rainfall-runoff factor, R, in the USLE with a runoff energy factor to 
produce the MUSLE, 


= O50 % nantes XK Suleiaciceeiin 


where S is sediment yield in tons, Q is runoff volume in acre-feet, 

q, is peak flow rate in ee /sec, and K LS C P are as defined in the 
USLE. The coefficient, 95, and exponent, 0.56, were determined by 

optimization, Williams, 1975. 


Runoff and suspended sediment yield data from four Reynolds Creek 
Experimental Watershed stations at Outlet, Macks Creek, Tollgate, and 
Reynolds Mountain, Figure 1, were analyzed using the MUSLE, Table 
4.5. First, data from individual events (runoff hydrographs and 
Sediment transport records) were separated into rainfall, snowmelt, 
and mixed, rainfall and snowmelt categories. Next, data for the 


Table 4.5.--Watersheds used in MUSLE analysis. 








Watershed Drainage Mean Average Number 
area elevation slope of events 
-acres- a——ftaa— ---%--- 
Reynolds Outlet 57,754 4900 a7 2242 
Macks Creek 7,846 4930 18 160 
Reynolds Tollgate 13,453 6100 26 330 
Reynolds Mountain 100 6830 14 371 
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following runoff and sediment yield components of each event were 
determined: 


1. Total peak streamflow rate in ft>/sec. 
2. Total runoff volume in acre-feet. 
3. Total sediment yield in tons. 


4. Storm peak streamflow rate in ft°/sec (total peak flow rate 
less the base flow rate). 


5. Storm runoff volume in acre-feet (total volume less the 
baseflow volume). 


6. Storm sediment yield in tons (total sediment yield less the 
sediment yield associated with base flow). 


Minimum, maximum, and mean values in each category are summarized in 
Table 4.6. Generally, events included a wide range of values with 
mixed rain and snowmelt contributing greatest peak flows, volumes, 
and sediment yields. The optimized best fit equations for each 
Watershed and category are listed in Table 4.7, assuming constant 
KLSCP values. Overall, storm runoff components showed a slightly 
better fit than total runoff components. Equations from rainfall 
events were most similar to the Williams' equation. Multiple 
correlation coefficients, r, ranged widely with Macks Creek data 
showing best results. 


Average yearly suspended sediment yields are shown in Table 4.8, to 
evaluate the relative importance of rainfall, snowmelt, and mixed 
rain and snowmelt in transporting sediment from each watershed. The 
Reynolds Outlet and Macks Creek watersheds, mean elevation about 4900 
feet, show nearly 90 percent sediment yield from mixed events. 
Reynolds Creek at Tollgate, mean elevation 6100 feet, shows 63 
percent sediment yield from mixed events, and 36 percent from 
snowmelt. Reynolds Mountain, mean elevation 6800 feet, shows 66 
percent sediment yield from snowmelt. 


Rangeland Soil Loss and Sediment Yield Predictions: 

The USLE has not been thoroughly validated on rangelands. 
Consequently, it should be used mainly to show relative effects of 
fires, grazing, seeding, land disturbance, and other land use changes 
on small areas where sediment deposition is not significant. Major 
questions on effects of snowmelt and frozen soil runoff, rocky soils, 
Slope length, extremely steep slopes, and rangeland vegetation are 
yet to be answered in accurately predicting soil losses. Generally, 
conventional application of the USLE appears to overpredict rangeland 
soil losses, compared with measured sediment yields. 


Results of predicting sediment yields using the PSIAC procedure show 
reasonable agreement with measured yields from sagebrush rangeland 
watersheds. The procedure is relatively simple, can be applied to 
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Table 4.6.--Summary of peak flow rates, runoff volumes, and sediment yields for the 
period of record, Reynolds Creek Experimental Watershed. 


C—O 





Watershed Category Range Peak streamflow Runoff volume Sediment vield 
and period and no. of Total Storm Total Storm Total Storm 
of record of events values ---ft>/sec---- --acre-feet--  # ----- tons----- 

Reynolds Rainfall Min. yl 28 8 5 20 19 

Outlet (11) Max. 821 815 200 Si, 1881 1879 

(1967-80) Mean 187 146 117 33 321 302 

Snowmelt Min. 39 11 47 8 72 4 

(91) Max. 225 126 449 112 eee 852 

Mean 5 i BP 48 178 39 139 92 

Mixed Min. 17 2 ahd al at Oeste 

(140) Max. 1663 1653 2689 2032 34713 33766 

Mean 198 127 284 108 1116 968 

Macks Rainfall Min. 3.9 3.6 2 1.4 1 0.9 

Creek (4:3) Max. 86 79 129 51 268 268 

(1968-80) | Mean 32 28 28 15 66 64 

Snowmelt Min. 3 2 3 1 Ore O24 

(58) Max. 71 50 123 TES 313 303 

Mean ou 13 25 9 ae 19 

Mixed Min. 3 1 4 1 Os 0.1 

(89) Max. SY27/ 304 486 357 5501 Syskone 

Mean 53 41 64 a2 288 274 

Reynolds Rainfall Min. abs Yor wait ae aby 8 5 5 

Tollgate (5) Max. 102 93 GES 30 46 44 

(1967-80) Mean 43 39 30 18 a2 20 

Snowmelt Min. 16 5 24 4 4 " 

(190) Max. 218 100 337, 82 766 684 

Mean 89 28 144 24 98 69 

Mixed Min. ) i 11 1 il 0.5 

(135) Max. 403 345 1012 616 7330 7225 

Mean 100 52 170 By} 256 230 

Reynolds Rainfall Min. O22 0.2 0.2 0.1 OL Og 

Mountain C2) Max. LoD 0.9 alae 0.3 00. -05 

(1969-80) Mean 0.9 0.6 Oa7 0.2 -03 -03 
Snowmelt Min. O.f 0.02 0.1 0.02 -002 -O01 

(259) Max. Soo 4.4 7.7 3.40 224 2a 
Mean Dee: 0-9 By ORF 0.24 0.18 
Mixed Min. On Orn 0.4 Oe 0.01 -001 
(110) Max. Cee Gre 93 4D 2.6 Vee 
Mean ALR Med Bisa, 0.83 0.29 e222 
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Table 4,7,--MUSLE optimization fitted coefficients and exponents, Reynolds 
Creek Experimental Watersheds. 


Watershed Category 
Reynolds Rainfall total 
ee Snowmelt total 
Mixed total 
Rainfall storm 
Snowmelt storm 
Mixed storm 
Macks Rainfall total 
Reeek Snowmelt total 
Mixed total 
Rainfall storm 
Snowmelt storm 
Mixed storm 
Reynolds Rainfall total 
Boles ce Snowmelt total ™ 
Mixed total 
Rainfall storm 
Snowmelt storm 
Mixed storm 
Reynolds Snowmelt total 
Mountain 


Mixed total 
Snowmelt storm 


Mixed storm 








Equation KLSCP og 
s = 0.01 (Qxq,)1"* 0.0166 0.79 
s = 0.02 (Qxqp)!*? 0.0166 0.59 
S$ = 1.02 (Qxqp) 9°97 0.0166 0.85 
s = 2.28 (Qxq,)1°™ 0.0166 0.99 
gspvo2 icgiqs) Mie 0.0166 0.67 
s = 266 (Qxq,)°-°8 0.0166 0.80 
s = 234 (Qxqp)9-4? 0.0183 0.71 
s = 0.22 (Qxq,) 130 0.0183 0.98 
s = 4.05 (Qxqp)9°% 0.0183 0.97 
s= 178 (Gaga) oo 0.0183 0.82 
s = 1.92 (Qxq,)1*?? 0.0183 0.96 
s = 26.5 ee a 0.0183 0.97 
s = 39.1 (Qxq,)°°?° 0.0171 0.84 
s = 0.05 (Qxqy)**7? 0.0171 0.73 
s = 0.02 (Qxq,) 1+? 0.0171 0.94 
s = 54.1 (Qxq,) 0.0171 0.79 
s = 7.64 (Qxqp)°*7" 0.0171 0.78 
s = 0.10 (Qxqp)2-7° 0.0171 0.91 
s = 1.67 (Qxqp) °°? 0.0221 0.63 
s = 3.48 (Qxqp)°°°7 00221 0.60 
s'= 7.42 (Qxa,)°" 0.0221 0.80 
s = 11.45(Qxqp)°*?? 0.0221 0.58 


Table 4.8.--Average yearly total suspended sediment yields by runoff category, 


Reynolds Creek Watersheds. 





Watershed Sediment yield, tons (percent of total) 
Rainfall Snowmelt Mixed rain-snow Total 
Reynolds Outlet 589 (4) 1,149(8) 12 ,014(88) UE US 
Macks Creek 143(6) 128(5) 2,331(89) 2,602 
Reynolds Tollyate 37(1) 1,545 (36) 2,662(63) 4,244 
Reynolds Mountain 0.03(0) 5.26(66) 2.0734) 7.96 
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large areas, and has factors to account for range treatments and 
management changes. It is most useful in predicting long-term annual 
yields and in identifying problem areas. 


The MUSLE has great potential for use in rangeland modeling, where 
data are available to accurately evaluate individual factors and 
interactions among processes associated with sediment movement are 
understood. Applications of the equation on an event basis encourage 
its increased use in modeling large watershed areas. 
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5. WATER QUALITY 


Personnel Involved 


G. R. Stephenson, Responsible for coordinating 

Geologist activities with cooperators. 
Designs collection network and 
responsible for project completion. 





KR, €. Rychert, Cooperator in microbiological 
Microbiologist aspects of rangeland streams as 
(Summer Employment) affected by grazing practices. 
A. P. Veigel, Responsible for collection of 
Seientific Aid Water samples and laboratory 
(U. of Idaho Cooperator ) analyses; assists in data 
processing. 
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WATER QUALITY OBSERVATIONS ~ BOISE FRONT 


Water quality baseline information continues to be developed from 
Samples collected at four sites on the Boise Front rest-rotation 
grazing system. The four sites are located immediately upstream from 
weirs, which record runoff. The sites are located such that the data 
Will reflect the grazing practices. Figure 5.1 gives the site 
location and the scheduled grazing activity for 1981. 


9/31/81 








jem emmennoaa| 
SCALE IN MILES 
| LEGENO 
nko te Past URE UCUHoARt —<—— FALL SHEEP DRIVE (:980) 
hts eke ——> SPRING SHEEP DRIVE (1981) 
t SHAE eaten wists) CATTLE GRAZING 
H HIGH PASTURES Nee Sea L4 - 4/15 TO 5/27/31 
1 * Romais) TOS IES, re ATTLE GRAZING 
. SS GC L 
auobaues L3 - 5/28 TO 6/25/31 
HS =" 6/7 25u LO ays St 





Figure 5.1.--Boise Front site locations and grazing activity, i951. 
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Table 5.1 gives the water quality data for the four sites for 


1979-1981 water years (Oct. 1 - Sept. 30). Total number of samples 
varies because of intermittent streamflow at several sites. 


Table 5.1.--Water quality characteristics, Boise Front Watershed sampling sites, 1979-81. 








No. of 
Parameters Units Samples = -=-=----~--— Max imum--—------- ---------Minimum------- --- -—----------Average-------- == 
(1981) 1979 1980 1981 1979 1980 1981 1979 1980 1981 








Os oS Ir aR 11S. ees Bk eS Rell 





pH units 18 8.35 7.95 8.00 6.67 7.20 7.10 7.73 7.61 7.54 
Conductivity umhos 15 216.00 207.62 252.88 11LE. 23 89.59 133.95 171.51 136.49 176.50 
Dissolved solids mg/l 15 149.97 143.93 174.48 54.53 62.46 92.43 114.05 95.0) 121.83 
Calcium mg/1 18 19.34 17.90 20.80 10.52 7.63 9.90 14.52 12.68 15.58 
Magnesium mg/l 18 2.95 2.70 3.90 2.00 1.66 2.15 2.35 2.14 3.03 
Sodiua mg/l 17 15.49 10.80 9.80 3.00 7.30 6.70 7.90 8.95 8.64 
Phosphorous mg/L 13 0.03 0.27 0.16 0.03 0.03 0.05 0.03 0.09 0.038 
Nitrate mg/L 6 0.03 1.16 0.40 0.03 0.64 0.28 0.03 0.84 0.35 
$io, mg/1 -- 36.00 40.00 -- 21.50 27.00 = 32.78 33.50 -- 
Sodium absorption ratio ratio — — — -- -- -- -- -- -- -- 
Suspended solids mg/1 5 167.60 124.00 20.00 3.00 3.00 0.00 21.28 33.89 8.60 
Total coliform cts/100 ml 16 980.00 38400.00 120000.00 0.00 300.00 600.00 234.00 8139.31 22891.25 
Fecal coliform ets/100 ml 17 675.00 580.00 3350.00 0.00 0.00 5.00 81.52 49.07 252555 
Fecal strep cts/100 ml 18 3320.00 3700.00 1360.00 0.00 8.00 4.00 431.81 473.86 298.72 
cov mg/L om 6.14 27.78 -- 0.00 0.00 -— 2.89 6.84 = 
BOD mg/1 1 3.40 1.90 5.50 0.00 0.00 5.50 1.05 1.75 Sac, 
DO ng/1 17 10.00 10.50 10.00 5.50 7.50 5.00 8.35 9.13 7.94 
Camp Creek 
pH units 9 8.20 8.00 7.80 6.92 7.40 7.20 7.68 7.72 7.58 
Conductivity umhos 6 166.80 200.21 334.60 142.80 129.20 167.00 153.78 161.24 210.56 
Dissolved solids mg/1 6 115.92 139.10 230.87 99.36 89.73 115.23 107.02 111.54 3465.29 
Calcium mg/L 9 11.23 18.00 27.40 12.23 10.80 12.80 11.23 2 epee as! 17.69 
Magnesium ng/1 9 2.64 3.40 8.90 2.64 2.42 1.60 2.64 2.86 4.29 
Sodium ng/1 9 9.76 10.90 11.30 9.76 7.90 7.40 9.76 8.94 anas 
Phosphorous mg/l 7 0.03 0.23 0.18 0.03 0.04 0.08 0.03 0.10 0.12 
Nitrace mg/1 1 0.03 -- 0.14 0.03 —- 0.41 0.03 -- 0.41 
Si0, mg/1 -— 32.00 45.00 -- 32.00 27.00 - 32.00 36.00 -- 
Sodium absorption ratio ratio 3 od -- 9.20 -- -- 3.00 -—- -- 6.07 
Suspended solids mg/L = 19.00 161.50 — 3.00 3.00 -- 11.50 53.46 -- 
Total coliform cts/100 ml 8 775.00 22000.00 170000 .00 12.00 20.00 2500.00 139.67 4946.80 63625.00 
Fecal coliform ects/100 ml 7 528.00 270.00 13000.00 0.00 0.00 4.00 110.33 45.72 4062.28 
Fecal strep cts/100 ml 9 550.00 2000.00 1800.00 4.00 6.00 8.00 124.33 248.08 751.11 
coD mg/1 = — 26.11 - -- 0.00 _ oe 8.78 -- 
BOD mg/1 _ == 1.50 -- -_ 0.50 = = 1.00 -- 
Do mg/1 8 10.00 10.50 10.00 9.50 7.00 5550 9.67 8.54 7-19 


zs Uppec Maynacd 


pH units 26 8.20 8.20 8.00 6.53 7.30 7.10 7.70 7.70 7.61 
Conductivity umhos 26 227.05 199.26 280.00 96.00 78.40 129.20 £59.32 153.33 139.25 
Dissolved solids me/L 26 158.63 167.81 193.20 66.65 54.43 89.15 110.77 103.25 123.68 
Calcium my/t 27 18.61 27.90 28.50 9.10 6.47 5.54 EZ on 14.61 16.73 
Magnesium mg/L oa 2.87 3.68 5.85 1.65 1.32 1.56 1.92 2567 3.046 
Sood fuera my /t 24 R00 12.80 12.10 2.90 6.25 6.00 3.65 JLt4 Roo 
Phosphorous mg/1 24 0.03 0.22 0.19 0.03 0.03 0.04 0.02 9.09 0.09 
Nitrate me/1 6 0.02 1.41 1 ee 0.02 0.66 0.25 0.02 1.03 0.56 
sia, mg/L — 39°55 37.00 —= 20.00 26.75 oer 29.78 30.42 == 
Sod tum absorption racio tatio -- =— -- -- -- - —— — -- -- 
Suspended solids mg/l 5 23.50 120.40 9.00 2.00 3.00 2.80 7.63 30.58 3.56 
Total coliform cts/100 ml 22 2300.00 40000 .00 268000.00 0.00 100.00 1500.00 382.13 9833.43 36950.00 
Fecal coliform cts/100 mL 26 1322.00 3820.00 20600 .00 0.00 0.00 9.00 230.13 310.27 1262.00 
Fecal strep cts/100 ml 26 3220.00 9300.00 7600.00 13.00 18.00 5.00 368.52 919.78 853.50 
cop me/1 5 13.10 25.19 8.380 0.00 0.00 1.70 4.80 7.38 4.96 
BUD mg/l Ll 2.50 2.50 4.80 0.00 0.00 4.80 0.55 0.96 4.80 
DO mg/L 28 10.00 10.50 12.00 7.00 5.00 5.00 8.24 8.50 8.13 


Highland Valley 


pH units 42 8.10 9.20 8.00 6.50 1h We 7.00 7.52 7.62 7.45 
Conductivity umhos 4L 275.63 232.00 310.00 118.30 118.80 111.78 753 149.90 173.06 
Disselved solids mg/1 41 192.51 161.46 213.90 BZ 511 75.62 77.13 119.38 107.30 119.41 
Calcium mg/L 43 22.37 16.50 19.10 10.79 8.20 5.39 17.43 13.30 12.55 
Magnesium mg/l 43 3.15 3.58 7.05 2.30 1.86 2.62 2.83 2.92 4.50 
Sudiunm mg/1 42 14.28 8.80 15.70 2.80 6.10 6.60 9.52 7.60 8.51 
Phosphorous mg/l 35 0.16 0.30 0.21 0.05 0.07 0.06 0.10 0.14 0.13 
Nitrate mg/L 14 0.24 6.33 9.16 : 0.24 2.82 Lore 0.24 4.37 4.08 
$10, mg/l -_ 42.21 42.00 -- 26.50 27.50 i 34.61 34.42 -- 
Sodium absorption ratio ratlo -- =< -- -- ~=— =< — peer a -- 
Suspenued solids me/1 15 210.70 621.80 178.00 2.00 14.50 4.00 23.76 LPT ..92 31 37 
Total colirorm ets/100 ml 39 31000.00 32000.00 490000 .00 0.00 100.00 1000.00 3646.75 8308.72 39132.82 
Fecal coliform cts/100 ml 40 23200.00 3500.00 6280.00 0.00 14.00 0.00 1295.06 240.66 704.75 
Fecal strep ets/100 mb 41 99000.00 9300.00 11800 .00 0.00 12.00 15.00 5655.69 1020.47 1333.78 
cop mg/L 6 39.57 92.70 12.70 2.00 0.58 2.40 9.61 14.42 7.10 
BOD 2 mg/l 9 6.40 5.00 9.20 0.00 0.00 4.70 1.65 1.67 7.49 
DO mg/l 43 9.50 10.50 10.50 6.00 6.00 5.50 7.76 8.66 7.85 
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When comparing the data from Table 5.1 for 1979-1981, the most 
Significant indicator of water quality variations continues to be the 
fecal coliform concentration. The fecal coliform indicator relates 
directly to the presence of warm blooded animals and reflects, 
mainly, the effects of the rest-rotation grazing management system. 
Other indicators, such as suspended solids and nutrients, show some 
variation between sites and between sampling dates, but are more 
frequently the result of runoff characteristics than the absence or 
presence of livestock. 


Table 5.2 gives the average fecal coliform concentration for the four 

Sites for the past three years. Each site reflects the water quality 

of the runoff from the individual grazing pastures within the 

| rest-rotation system. During the years cattle graze the individual 
pastures, the average fecal coliform concentration is usually the 
highest. In 1980, pasture L2, the Lower Maynard site, was grazed, 
but the stream was dry when the cattle were turned in, so no Samples 





| were taken. The cattle obtained water from spring flow developments. 
| In 1981, pasture H2, Upper Maynard, was not scheduled for grazing, 
but had a very high fecal coliform average. The reason for this was 

that, because of low stream flow, gates between pastures were left 
open so cattle could move between pastures H3 and H2. Consequently, 
| pasture H2 was grazed in 1981 even though the grazing schedule called 
| for the pasture to be rested. 


Table 5.2.--Average fecal coliform concentrations. 


| 





| 1978 1979 1980 1981 
| Lower Maynard (12) 225; ; 82 “49 °** 252 
| Camp Creek (L3) 38 110 45 4062 | 
| Upper Maynard (H2) 21g 230 310 1262 
| Highland Valley (H3) 193 izes 240 704 


*Framed values indicate field grazed that year. 





**Cattle not turned in until after October 1, 1980. Streams were dry; 
no samples. 
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Figure 5.2 shows graphically the variations in fecal coliform 
concentration throughout the water year at all four sites on the 
Boise Front. Evidence of the large deer herd wintering in this area 
is recorded in periodic increases in fecal coliform concentration at 
all sites through the winter and spring months. The Highland Valley 
Site (H3) shows the largest influence of deer, especially in April 
and May, as the herd moves through tnat area to higher country for 
the summer. 
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Figure 5.2.--Fecal ‘coliform concentration, Boise Front sampling sites. 
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The 1980 fall sheep drive passed through pastures H3, L3, and H2 in 
November. Evidence of their passage, recorded in elevated fecal 
coliform concentration at the stream site, is recorded only at the 
Highland Valley and Upper Maynard sites. The Camp Creek site, 
located in pasture L3, was dry. The 1981 spring sheep drive occurred 
during May, and passed through L2, L3, H3, and H2, respectively. 
Elevated fecal coliform counts are recorded at all sampling sites in 
these pastures at this time. Cattle entered low pasture 3 in late 
May. The effect of their presence was recorded immediately in 
elevated fecal coliform counts at the Camp Creek site as they 
frequented the riparian habitat for water, shade, and more palatable 
grass. The cattle moved out of pasture L3 into pasture H3 in late 
June, but fecal coliform concentration remained high until the stream 
above the Camp Creek site dried up in early July. 


Even though cattle were scheduled to be in pasture H3 by early July, 
the streams were running low and gates into L3, L2, and H2 were 
opened so the cattle could utilize these pastures. As the weather 
became hotter and drier, all the streams dried up, except above the 
Highland Valley site, high pasture 3. As seen in ri gure 85 .2.,"'the 
elevated fecal coliform concentration at this site remained high 
until late summer as the cattle moved into this area for better 
grass, Water, and cooler air temperatures. 


Model evaluation and parameter testing: 


Water quality model guidelines were evaluated this year. Since 
bacterial parameters appear to be the most sensitive indicator of 
nonpoint source pollution in rangeland streams, attempts are being 
made to develop a user oriented microbiology predictive method to 
Simulate bacterial variations in these environments. During this 
first attempt, the following independent variables were selected 
along with fecal coliform concentration as the dependent variable: 
runoff, dissolved oxygen, water temperature, average slope, aspect, 
drainage density, vegetative characteristics, channel roughness, and 
grazing practice (number of cattle). Relationships between land form 
factors, such as slope steepness and drainage density, and bacterial 
factors were found to exist in most grazing allotments. The greater 
the drainage density (stream length per unit area), the lower the 
overall bacterial concentration; and the steeper the slopes, the 
higher the overall bacterial concentration. Both factors influence 
the distribution of livestock and their proximity to streams. 


Regression techniques were used to test for parameter relationships. 
The initial test, using dissolved oxygen, temperature, and streamflow 
as independent variables, and fecal coliform concentration as the 
dependent variable proved unsuccessful, as correlations were not 
Significant. However, individual parameter correlations with fecal 
coliforms, using linear regressions, were significant. Suspended 
sediment, r@ = -/9, WaS correlated with fecal coliforms, as was 
channel roughness, r* = .96. The latter parameter had a negative 


Slope. Testing of these parameters will continue for model 
development. 


Fats 


E. coli occurrence and survival: 

Several detailed studies were completed this year to determine the 
survival of E. coli in different environments. In order to test 
previous findings, which demonstrated that stream sediments can serve 
as a reservoir of E. coli, a detailed sampling network was carried 
out during January through October. The mean E. coli concentrations ~ 
for the sediment to water ratio was 13.5, confirming previous 
findings. E. coli survival in the stream sediments is often carried 
over from one grazing season to another, and frequently experiences 
regrowth. Ee coli in the overlying stream water, however, seldom 
remain viable in the water for any length of time. The stream 
environment varies considerably with roughness, temperature, and 
predators, etc.; whereas, the sediment environment is quite stable. 
Disturbance of the sediment causes resuspension of the bacterial 
indicators in the overlying waters. 








In order to determine seasonal die-off, regrowth, and movement of Ee 
coli (fecal coliforms) in fresh fecal material in rangeland 
environments, eight paired sites were selected with varying slopes, 
aspects, vegetation, and soil characteristics. Weekly samples were 
taken from cow pies freshly deposited in May and in October. Results 
show that E. coli populations die off rapidly during the summer 
months, because of dessication of the fecal material. Less than 0.1 
percent of the original population remained after two months. Weekly 
Samples from fresh October deposits, subject to freezing and thawing, 
show a die-off rate of one order of magnitude less, with a low level 
of regrowth occuring on south facing slopes during the following 
Spring. This is followed by complete summer die-off. Movement of E. 
coli from the cow pies into the soil following rain storms was <n 
negligible, even on the steeper slopes. Slope, aspect, soil type, 
and vegetative composition show no consistent relationships with 
die-off or regrowth. Final conclusions of these experiments are: 

(1) fecal material from livestock, as sources of fecal bacteria in 
rangeland streams, must either be deposited in the stream or 
transported to the streams by animals or immediately adjacent runoff; 
and (2) E. coli (fecal coliforms) die off rapidly during hot dry 
weather, but may remain viable during cool, moist weather, and can 
experience regrowth the following spring. 











In the continuing investigation of the cause of atypical E. coli in 
rangeland streams, the effects of antibiotics in cattle feed were 
investigated. Preliminary results indicate that 13 percent of the 
atypical E. coli isolates were resistant to two or more antibiotics,; 
whereas, 72 percent of the typical E. coli were resistant. This 
indicates that a higher percentage of the fecal coliforms in these 
rangeland streams are multiply-resistant to antibiotics, and could 
represent a significant health hazard not previously recognized. 
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